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Introduction 



This text is intended as an introduction to magnetohydrodynamics in astrophysics, 
emphasizing a fast path to the elements essential for physical understanding. It assumes 
experience with concepts from fluid mechanics: the fluid equation of motion and the La- 
grangian and Eulerian descriptions of fluid flow. In addition, the basics of vector calculus 
and elementary special relativity are needed. Not much knowledge of electromagnetic 
theory is required. In fact, since MHD is much closer in spirit to fluid mechanics than 
to electromagnetism, an important part of the learning curve is to overcome intuitions 
based on the vacuum electrodynamics of one's high school days. 

The first chapter (only 36 pp) is meant as a practical introduction including exer- 
cises. This is the 'essential' part. The exercises are important as illustrations of the 
points made in the text (especially the less intuitive ones). Almost all are mathemati- 
cally unchallenging. The supplement in chapter 2 contains further explanations, more 
specialized topics and connections to the occasional topic somewhat outside MHD. 

The basic astrophysical applications of MHD were developed from the 1950s through 
the 1980s. The experience with MHD that developed in this way has tended to remain 
confined to somewhat specialized communities in stellar astrophysics. The advent of 
powerful tools for numerical simulation of the MHD equations has enabled application 
to much more realistic astrophysical problems than could be addressed before. The 
availability of these tools has also made magnetic fields attractive to a wider community. 
In the course of this development the interpretation of numerical results, as well as of 
magnetohydrodynamics in general has suffered somewhat. 

This text aims to show how MHD can be used more convincingly when armed with 
a good grasp of its intrinsic power and peculiarities, as distinct from those of vacuum 
electrodynamics or plasma physics. The emphasis is on physical understanding by the 
visualization of MHD processes, as opposed to more formal approaches. This allows one 
to formulate better questions more quickly, and to interpret computational results more 
meaningfully. For more comprehensive introductions to astrophysical MHD, see Parker 
(1979), Kulsrud (2005) and Mestel (2012). 

In keeping with common astrophysical practice Gaussian units are used. 



1 



1. Essentials 



Magnetohydrodynamics describes electrically conducting fluidq^ in which a magnetic 
field is present. A high electrical conductivity is ubiquitous in astrophysical objects. 
Many astrophysical phenomena are influenced by the presence of magnetic fields, or 
even explainable only in terms of magnetohydrodynamic processes. The atmospheres of 
planets are an exception. Much of the intuition we have for ordinary earth-based fluids 
is relevant for MHD as well, but more theoretical experience is needed to develop a feel 
for what is specific to MHD. The aim of this text is to provide the means to develop 
this intuition, illustrated with a number of simple examples and warnings for popular 
pitfalls. 



1.1. Equations 

1.1.1. The MHD approximation 

The equations of magnetohydrodynamics are a reduction of the equations of fiuid me- 
chanics coupled with Maxwell's equations. Compared with plasma physics in general, 
MHD is a strongly reduced theory. Of the formal apparatus of vacuum electrodynamics, 
with its two EM vector fields, currents and charge densities, MHD can be described 
with only a single additional vector: the magnetic field. The 'MHD approximation' that 
makes this possible involves some assumptions: 

1. The fluid approximation: local thermodynamic quantities can be meaningfully de- 
flned in the plasma, and variations in these quantities are slow compared with the time 
scale of the microscopic processes in the plasma. 

2. In the plasma there is a local, instantaneous relation between electric fleld and cur- 
rent density (an 'Ohm's law'). 

3. The plasma is electrically neutral. 

This statement of the approximation is somewhat imprecise. I return to it in some of 
the supplementary sections of the text (chapter [2]) . The first of the assumptions involves 
the same approximation as used in deriving the equations of fiuid mechanics and ther- 
modynamics from statistical physics. It is assumed that a sufficiently large number of 
particles is present so that local fluid properties, such as pressure, density and velocity 
can be defined. It is sufficient that particle distribution functions can be defined prop- 
erly on the length and time scales of interest. It is, for example, not necessary that the 
distribution functions are thermal, isotropic, or in equilibrium, as long as they change 
sufficiently slowly. 

The second assumption can be relaxed, and the third is closely related to the second 



(cf. sect. 2.6). For the moment we consider them as separate. The second assumes that 
whatever plasma physics takes place on small scales, these average out to an instanta- 
neous, mean relation (not necessarily linear) between the local electric field and current 
density, on the length and time scales of interest. The third assumption of electrical 



^ In astrophysics 'fluid' is used as a generic term for a gas, liquid or plasma 
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neutrality is satisfied in most astrophysical environments, but it excludes near- vacuum 



conditions such as the magnetosphere of a pulsar (sect. 2.10). 

Electrical conduction, in most cases, is due to the (partial) ionization of a plasma. 
The degree of ionization needed for 2. to hold is generally not large in astrophysics. The 
approximation that the density of charge carriers is large enough that the fluid has very 
little electrical resistance: the assumption of perfect conductivity, is usually a good first 
step. Exceptions are, for example, pulsar magnetospheres, dense molecular clouds or the 
atmospheres of planets. 



1.1.2. Ideal MHD 

Consider the MHD of a perfectly conducting fluid, i.e. in the limit of zero electrical 
resistance. This limit is also called ideal MHD. Modifications when the conductivity is 
finite are discussed in sections 11.101 and 12.91 

The electric field in a perfect conductor is trivial: it vanishes, since the electric current 
would become arbitrarily large if it did not. However, the fluid we are considering is 
generally in motion. Because of the magnetic field present, the electric field vanishes 
only in a frame of reference moving with the fiow; in any other frame there is an electric 
field to be accounted for. 

Assume the fluid to move with velocity v(r) relative to the observer. Let E' and 
B' be the electric and magnetic field strengths measured in an instantaneous inertial 
frame where the fluid is at rest (locally at the point r, at time t). We call this the 
comoving frame or fluid frame. They are related to the fields E, B measured in the 
observer's frame by a Lorentz transformation (e.g. Jackson E&M Chapter 11.10). Let 
= V • E/u be the component of E parallel to the fiow, E_|_ = E — vE • v/v"^ the 
perpendicular component of E, and similar for B. The transformation is theiJ^ 



E'^ = 7(Ex + vxB/c), 
B'w = , 

B'^ = 7(Bx-vxE/c), (1.1) 

where 7 is the Lorentz factor 7 = (1 — v"^ /c'^)~^^'^ , and c the speed of light. By the 
assumption of infinite conductivity, E' = 0. The electric field measured by the observer 
then follows from (1.1) as: 

E = -vxB/c. (1.2) 

Actually measuring this electric field would require some planning, since it can be 
observed only in an environment that is not itself conducting (or else the electric field 
would be shunted out there as well). The experimenter's electroscope would have to be 
kept in an insulating environment separated from the plasma. In astrophysical appli- 
cation, this means that electric fields in MHD become physically significant only at a 
boundary with a non-conducting medium. The electric fields associated with differential 
flow speeds within the fluid are of no consequence, since the fluid elements themselves 
do not sense them. 

A useful assumption is that the magnetic permeability and dielectrical properties of 
the fluid can be ignored, a very good approximation for many astrophysical applications. 



^ This transformation is reproduced incorrectly in some texts on MHD. 
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This allows the distinction between magnetic field strength and magnetic induction, 
and between electric field and displacement to be ignored. This is not an essential 
assumption. Maxwell's equations are then 

47rj + d^/dt = cV X B, (1.3) 

d'B/dt = -cV X E, (1.4) 

V-E = 47rcj, (1.5) 

V-B = 0, (1.6) 

where j is the electrical current density and a the charge density. Taking the divergence 

onservation of charge: 

+ V-j = 0. (1.7) 



of Maxwell's equation (1.3) yields the conservation of charge: 

da 



dt 

1.1.3. The induction equation 



Using expression (1.2) for the electric field, the induction equation (1.4) becomes 

— = Vx(vxB). (1.8) 

This is known as the induction equation of ideal MHD, or MHD induction equation 
for short. It describes how the magnetic field in a perfectly conducting fluid changes 



with time under the influence of a velocity field v (see section 1.10 for its extension to 
cases when conductivity is finite). 

By the MHD approximation made, the original electromagnetic induction equation 
has changed flavor drastically. From something describing the generation of voltages by 
changing magnetic fields in coils, it has become an evolution equation for a magnetic 
field embedded in a fluid flow. For historical reasons, it has retained the name induction 
equation even though it actually is better understood as something new altogether. The 



divergence of (1.8) yields 

|v.B = 0^ (1.9) 

The MHD induction equation thus incorporates the condition V • B = 0. It need not 
be considered explicitly in solving the MHD equations, except as a requirement to be 
satisfied by the initial conditions. 

1.1.4. Geometrical meaning of V - B = 

The field lines of a divergence-free vector such as B (also called a solenoidal vector) 
'have no ends'. While electric field lines can be seen as starting and ending on charges, 
a magnetic field line (a path in space with tangent vector everywhere parallel to the 
magnetic field vector) wanders around without meeting such monopolar singularities. 
As a consequence, regions of reduced field strength cannot be local: magnetic field lines 
must 'pass around' them. This is illustrated in Fig. |l.l[ Unlike scalars like temperature or 
density, the magnetic field strength cannot be reduced or increased only locally. A change 
in field strength must be accommodated by changes in field line shape and strength in 



its surroundings (for an example and consequences see problem 3.20). 
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Figure 1.1.: Field lines near a field-free inclusion (dashed contour). The configuration on the 
left contains field lines ending on the surface of the inclusion. V • B = excludes 
such configurations. Regions of reduced field strength distort the field in their 
surroundings (middle). Right: pill box (red) to derive continuity of from div 
B = 0. 



A bit more formally, let S be a surface in the field configuration, at an arbitrary angle 
to the field lines, and let n be a unit vector normal to S. Apply Gauss's theorem to 
div B = in a thin box of thickness e — )• oriented parallel to S. The integral over the 
volume being equal to the integral of the normal component over the box then yields 
that the component of B normal to the surface, 

Bn = Bn (1.10) 

is continuous across any surface S. Hence also across the dashed surface in the left panel 
of Fig.pl 



1.1.5. Electrical current 



Up to this point, the derivation is still valid for arbitrary fluid velocities. In particular 



the induction equation (1.8) is valid relativistically, i.e. at all velocities v < c (though 
only in ideal MHD, not with finite resistivity). We now specialize to the nonrelativistic 
limit V <^ c. Quantities of first order in v/c have to be kept in taking the limit, since 



(1.1) 



the electric field is of this order, but higher orders are omitted. Substituting (1.2) into 
one finds that 

B[l + 0{v'/c% (1.11) 



B' 



i.e. the magnetic field strength does not depend on the frame of reference, in the non- 



relativistic limit. Substituting (1.2) into (1.3) yields: 



47rj 



d_ 
di 



V X B)/c = cV X B. 



(1.12) 



The second term on the left hand side, the displacement current, can be ignored if f ^ c. 
To see this, choose a length scale L that is of interest for the phenomenon to be studied. 
Then |V x B| is of the order B/L. Let ^ be a typical value of the fluid velocities 
relevant for the problem; the typical time scales of interest are then or order r = L/V. 
An upper limit to the displacement current for such length and time scales is thus of 
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order |v x B/t|/c ~ BiV / L){y / c)^ which vanishes to second order va. v/c compared to 
the right hand side. Thus (1.12) reduces tcj^ 



47r 



Vx B. 



Taking the divergence: 



0. 



(1.13) 



(1.14) 



Equation (1.14) shows that in MHD currents have no sources or sinks. As a conse- 



quence it is not necessary to worry 'where the currents close' in any particular solution 
of the MHD equations. The equations are automatically consistent with charge conser- 
vatior[^ It is not even necessary that the field computed is an accurate solution of the 
equations of motion and induction. As long as the MHD approximation holds and the 
field is physically realizable, i.e. V - B = 0, the current is just the curl of B (eq. 1.13), 
and its divergence consequently vanishes. 

Of course, this simplification only holds as long as the MHD approximation itself is 
valid. But whether that is the case or not is a different question altogether: it depends 
on things like the microscopic processes determining the conductivity of the plasma, not 
on global properties like the topology of the currents. 



1.1.6. Charge density 



With the equation for charge conservation (1.7), eq. (1.14) yields 



5a 



0. 



(1.15) 



We conclude that it is sufficient to specify a = in the initial state to guarantee that 
charges will remain absent, consistent with our assumption of a charge-neutral plasma. 
Note, however, that we have derived this only in the non-relativistic limit. The charge 



density needs closer attention in relativistic MHD, see section 2.6 



Eq. (1.15) only shows that a charge density cannot change in MHD, and one might 



ask what happens when a charge density is present in the initial conditions. In practice, 
such a charge density cannot not last very long. Due to the electrical conductivity of the 
plasma assumed in MHD, charge densities are quickly neutralized. They appear only at 
the boundaries of the volume where MHD holds. See 12.51 and 12.71 



1.1.7. Lorentz force, equation of motion 



With relation ( 1.13 ) between field strength and current density, valid in the non-relativistic 



limit, the Lorentz force acting per unit volume on the fiuid carrying the current is 



1. 



-j xB 



1 



(Vx B)xB. 



(1.16) 



See 



Perfect conductivity has been assumed here, but the result also applies at finite conductivity, 
problem 13.31 

The fact V • j = is stated colloquially as 'in MHD currents always close'. The phrase stems from 
the observation that lines of a solenoidal (divergence-free) vector in two dimensions have two choices: 
either they extend to infinity in both directions, or they form closed loops. In three dimensions it 
is more complicated: the lines of a solenoidal field enclosed in a finite volume are generally ergodic. 
A field line can loop around the surface of a torus, for example, never to return to the same point 
but instead filling the entire 2-dimensional surface. It is more accurate to say that since currents are 
automatically divergence-free in MHD, the closing of currents is not an issue. 
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This looks very different from the Lorentz force as explained in wikipedia. From a force 
acting on a charged particle orbiting in a magnetic field, it has become the force per 
unit volume exerted by a magnetic field on an electrically neutral, but conducting fluid. 

In many astrophysical applications viscosity can be ignored; we restrict attention here 
to such inviscid flow, since extension to a viscous fluid can be done in the same way as 
in ordinary fluid mechanics. Gravity is often important as an external force, however. 
Per unit volume, it is 

Fg = pg = -pV0, (1.17) 

where g is the acceleration of gravity, (p its potential and p the mass per unit volume. 
If p is the gas pressure, the equation of motion thus becomes 

dv 1 
P^ = -Vp+— (VxB)xB + pg, (1.18) 

where d/dt is the total or Lagrangian time-derivative, 

d/dt = d/dt + vV . (1.19) 

The Lorentz force in MHD is quadratic in B and does not depend on its sign. The 
induction equation (1.8): 

c?B 

— = V X (v X B) (1.20) 

is also invariant under a change of sign of B. The ideal MHD equations are therefore 
invariant under a change of sign of B: the fluid 'does not sense the sign of the magnetic 
fieldj^ Electrical forces do not appear in the equation of motion since charge densities 



are negligible in the non-relativistic limit (sect. 2.6). 

The remaining equations of fluid mechanics are as usual. In particular the continuity 
equation, expressing the conservation of mass: 

|^ + V-(pv) = 0, (1.21) 

or 

^ + pV-v = 0. (1.22) 

In addition to this an equation of state is needed: a relation p{p, T) between pressure, 
density, and temperature T. Finally an energy equation is needed if sources or sinks of 
thermal energy are present. It can be regarded as the equation determining the variation 
in time of temperature (or another convenient thermodynamic function of p and p). 



The equation of motion (1.18) and the induction equation (1.20) together determine 
the two vectors B and v. Compared with ordinary fluid mechanics, there is a new 
vector field, B. There is an additional equation for the evolution of this field: the MHD 
induction equation, and an additional force appears in the equation of motion. 

These equations can be solved without reference to the other quantities appearing 
in Maxwell's equations. This reduction vastly simplifies the understanding of magnetic 
fields in astrophysics. The price to be paid is that one has to give up most of the intuitive 
notions acquired from classical examples of electromagnetism, because MHD does not 
behave like EM anymore. It is a fluid theory, close in spirit to ordinary fluid mechanics 
and to the theory of elasticity. 



In non-ideal MHD this can be different, in particular when Hall drift is important (sect. 2.91 



7 



1.1.8. The status of currents in MHD 



Suppose V and B have been obtained as a solution of the equations of motion and 



induction (1.18, 1.20), for a particular problem. Then Ampere's law (1.13) can be used 
to calculate the current density at any point in the solution by taking the curl of the 
magnetic field. 

The nature of Ampere's law has changed: from an equation for the magnetic field 
produced by a current distribution, as in vacuum electrodynamics, it has been demoted 
to the status of an operator for evaluating a secondary quantity, the current. 

As will become apparent from the examples further on in the text, the secondary 
nature of currents in MHD is not just a mathematical curiosity. Currents are also rarely 
useful for physical understanding in MHD. They appear and disappear as the magnetic 
field geometry changes in the course of its interaction with the fluid flow. An example 



illustrating the transient nature of currents in MHD is given in problem 3.1 Regarding 
the currents as the source of the magnetic field, as is standard practice in laboratory 
electrodynamics and plasma physics, is meaningless in MHD. 

When familiarizing oneself with MHD one must set aside intuitions based on batteries, 
current wires, and induction coils. Thinking in terms of currents as the sources of B 
leads astray; 'there are no batteries in MHD'. (For the origin of currents in the absence 



of batteries see 2.8). Another source of confusion is that currents are not tied to the 
fiuid in the way household and laboratory currents are linked to copper wires. A popular 
mistake is to think of currents as entities that are carried around with the fluid. For 
the currents there is no equation like the continuity equation or the induction equation, 
however. They are not conserved in displacements of the fluid. 

1.1.9. Consistency of the MHD approximation 

In arriving at the MHD equations, we have so far accounted for 3 of the 4 of Maxwell's 



equations. The last one, eq. (1.5), is not needed anymore. It has been bypassed by 



the fact that the electric fleld in MHD follows directly from a frame transformation, eq. 



(1.2). Nevertheless, it is useful to check that the procedure followed has not introduced 
an inconsistency with Maxwell's equations, especially at relativistic velocities. This is 
done in section [2^ 



1.2. The motion of field lines 



In vacuum electrodynamics, fleld lines do not 'move' since they do not have identity that 
can be traced from one moment to the next. In ideal MHD they become traceable, as 
if they had an individual identity, because of their tight coupling to the fluid. (See 2.2 
about this coupling at the microscopic level.) 



This coupling is described by the induction equation (1.20). It does for the magnetic 



fleld (the flux density) what the continuity equation does for the mass density, but there 
are important differences because of the divergence-free vector nature of the fleld. To 



explore these differences, write the continuity equation (1.21) as 



dp 
di 



Vp. 



(1.23) 



The flrst term describes how the mass density p varies in time at some point in space 
as fluid contracts or expands. The second term, called the advection of the density p, 
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describes the change at a point in space as fluid of varying density passes by it. The 
induction equation can be written by expanding its right hand side, using the standard 



vector identities (sect. 4.1): 



dB 

Ik 



-B V • V - (v V)B + (B • V)v, 



(1.24) 



where V • B = has been used. In this form it has a tempting similarity to the continuity 
equation (1.23). The first term looks hke it describes the effect of compression and 



expansion of the fluid, hke the first term in (1.23) does for the mass density. The second 



term similarly suggests the effect of advection. There is, however, a third term, and as 
a consequence neither the effects of compression nor the advection of a magnetic field 
are properly described by the first and second terms alone. 

A form that is sometimes useful is obtained by combining the equations of continuity 
and induction: 

' ^ = ^ (.25, 



dt \ p / \ p 

This describes how the ratio of magnetic fiux to mass density changes when the fluid 



velocity varies along a fleld line (see problems 3.4 3.6) 




Figure 1.2.: A closed loop £ of fluid elements carried by the flow, with field lines passing 
through it. 



1.2.1. Magnetic flux 

While the induction equation does not work for B in the same way as the continuity 
equation does for the gas density, there is in fact a conservation of something playing a 
similar role, namely magnetic fiux. First we need to define magnetic fiux in this context. 

Consider a closed loop C of infinitesimal fiuid elements (Fig. 1.2). It moves with the 
fluid, changing its length and shape. The magnetic fiux of this loop is now defined as 
the 'number of fleld lines passing through' it. A bit more formally, let 5 be a surface 
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bounded by the loop. There are many such surfaces, it does not matter which one we 
take (problem 3.7). Then define the magnetic flux of the loop as 



$(£) = / B-dS, (1.26) 



s 

where dS = ndS*, with dS an element of the surface S and n the normal to S. In a 
perfectly conducting fluid the value of $ is then a property of the loop, constant in time, 
for any loop moving with the flow (Alfven's theorem): 

d^> 

-=0. (1.27) 



(slightly intuitively) in 2.1 



The equivalence of eq. (1.27) with the ideal MHD induction equation (1.20) is derived 



The flux ^{S) of the loop also defines a flux bundle: on account of V • B = the 
field lines enclosed by the loop can be extended in both directions away from the loop, 
tracing out a volume in space filled with magnetic field lines. Like the loop £ itself, this 
bundle of flux <I> moves with the flow, as if it had a physical identity. By dividing the 
loop into inflnitesimally small sub-loops, we can think of the fiux bundle as consisting of 
'single fiux lines', and go on to say that the induction equation describes how such flux 
lines move with the flow. This is also described by saying that held lines are 'frozen in' 
the fluid. Each of these fleld lines can be labeled in a time-independent way, for example 
by labeling the fluid elements on the surface S at some point in time to- 

If we define a 'fiuid element' intuitively as a microscopic volume carrying a fixed 
amount of mass (in the absence of diffusion of particles), a fiux line is a macroscopic 
string of such elements. It carries a fixed (infinitesimal) amount of magnetic fiux (in the 
absence of magnetic diffusion) through a cross section that varies along its length and 
in time. 



When the conductivity is finite, (1.27) does not hold. Field lines can 'diffuse across' 



the fiuid, and they cannot always be labeled in a time-independent way anymore. The 



induction equation then has an additional term describing such diffusion (sect. 1.10). 
Astrophysical conditions are, with some exceptions, close enough to perfect conductivity 
that intuition based on ideal MHD is applicable as a first step in most cases. The opposite 
limit of low conductivity is more amenable to applied mathematical analysis but rarely 
relevant in astrophysics. 

1.2.2. Field amplification by fluid flows 

In the absence of diffusion, magnetic fields embedded in a fiuid have a strong memory 
effect: their instantaneous strength and shape reflects the entire history of the fluid 
motions acting on it. One consequence of this is that fluid motions tend to rapidly 



increase the strength of an initially present magnetic fleld. This is illustrated in Fig. 1.3 



The fluid is assumed to have a constant density, i.e. it is incompressible. By the continuity 



equation (1.23) the velocity field is then divergence-free, V - v = 0. Assume that there 
is an initially uniform field in the z-direction, B = i?z. 

The field is deformed by fiuid motions which we assume to be confined to a volume of 
constant size. The velocities vanish at the boundary of this sphere. The field lines can 
then be labeled by their positions at the boundary of the sphere, for example a field line 
entering at xi, —z\ and exiting at xi, z\ can be numbered 1. Define the length L of the 
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Figure 1.3.: Field amplification by a complex fluid flow 



field line as its path length from entry point to exit; the initial length of line number 1 
is L = 2z\. 

Since the field lines are initially straight, any flow in the sphere will increase their 
length L. Let b be the average distance of field line number 1 from its nearest neighbor. 
As long as the field is frozen in, the mass enclosed between these field lines is constant, 
and on account of the constant density of the fluid, its volume Lb is also constanfj^ 



Hence b must decrease with increasing L as 5 ~ 1/L. By flux conservation (1.26), the 
field strength then increases as L. 



A bit more formally, consider the induction equation in the form (1.24), and by taking 
the second term on the right to the left hand side, write it as 

= _B V • v + (B -Vlv. (1.28) 

at 

It then describes the rate of change of B in a frame comoving with the fluid. Like the 
mass density, the field strength can change by compression or expansion of the volume 
(first term). This change is modified, however, by the second term which is intrinsically 
magnetic. Under the assumed incompressibility the first term vanishes, and the induction 
equation reduces to 

dB 

— = (B-V)v. (p = cst.) (1.29) 

To visualize what this means, consider two simple examples. In the first example we 
consider the effect of amplification of a magnetic field by a flow which converges on the 
field lines and diverges along them. Assume an initially uniform magnetic field in the 
x-direction (Fig.|1.4|) in a stationary, divergence-free velocity field: 



B(t = 0) = -Bqx, Vx = ax^ Vz = —az, Vy = 0. (1.30) 



At t=0, eq. (1.29) yields: 

dB, dB 



-,=o = ^U=o = i?oax. (1.31) 



^' This argument assumes a 2-dimensional flow, as suggested by Fig. 1.3 Except at special locations in 
a flow, it also holds in 3 dimensions. To see this requires a bit of visualization of a magnetic field in 
a shear flow. 
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Figure 1.4.: Field amplification by a converging-diverging flow. 
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Figure 1.5.: Field amplification by a shear flow perpendicular to B 



The first equality holds since at i = the field is uniform, so that (v •V)B = 0. Eq. 



( 1.31 ) shows that the field stays in the x-direction and remains uniform. The assumptions 



made at i = therefore continue to hold, and (1.31) remains valid at arbitrary time t 
and can be integrated with the result: 



Bq e"* X. 



:i.32) 



Stretching by a flow like ( 1.30 ) thus increases B exponentially for as long as the stationary 



velocity field is present. This velocity field is rather artificial, however. If the flow instead 
results from a force that stretches the left and right sides apart at a constant velocity, 
the amplification is only linear in time (exercise: verify this for yourself, using mass 
conservation and taking density constant). 

When the flow v is perpendicular to B everywhere, the term (B •V)v describes an 



effect that is conceptually rather different (see Fig. 1.5). Take the magnetic field initially 
uniform in the z-direction, the flow in the x-direction, with a piece-wise linear shear in 



z: 



and 



B{t = 0) = Bqz, 



(t = 0), 



(1.33) 



-vo (z < -1), 



VqZ (-1 < Z <1) 



iz>l] 



(1.34) 
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The induction equation (1.20) then yields (analogous to problem 



3.1): 



B,{t) = Bo= est, (1.35) 

= Bovot (-K z < 1), B^ = (z < -1, z > 1). (1.36) 

In the shear zone — 1 < z < 1 the x-component of B grows linearly with time. One could 
call this kind of field amplification by a flow perpendicular to it 'shear amplification'. 

Shear zones of arbitrary horizontal extent like the last example are somewhat artificial, 
but the effect takes place in essentially the same way in a rotating shear flow (such as 
a differentially rotating star for example). In this case, the field amplification can be 
described as a process of 'winding up' of field lines. A simple case is described by problem 

E33 

In these examples the flow field was assumed to be given, with the magnetic field 
responding passively following the induction equation. As the field strength increases, 
magnetic forces will eventually become important, and determine the further evolution 
of the field in a much more complex manner. Evolution of a magnetic field as described 
above, the so-called kinematic approximation, can hold only for a limited period of time 



in an initially sufficiently weak field (see also sect. 1.4) 



1.3. Magnetic force and magnetic stress 
1.3.1. Magnetic pressure and curvature force 

The Lorentz force is perpendicular to B. Along the magnetic field, the fluid motion is 
therefore subject only to the normal hydrodynamic forces. This makes the mechanics of 
a magnetized fluid extremely anisotropic. 



To get a better feel for magnetic forces, one can write the Lorentz force (1.16) in 



alternative forms. Using the vector identities (sect. 4.1) 



Fl = -^(V X B)xB = -^VB^ ^ J_(B.V)B. (1.37) 
47r ovr 47r 

The first term on the right is the gradient of what is called the magnetic pressure B'^ /Sir. 
The second term describes a force due to the variation of magnetic field strength in the 
direction of the field. It is often called the magnetic curvature force. 

These names are a bit misleading. Even in a field in which the Lorentz force vanishes. 



the magnetic pressure gradient (1.37) is generally nonzero, while the 'curvature' term 
can also be present when the field lines are straight. To show the role of curvature of 
the field lines more accurately, write the magnetic field as 

B = Bs, (1.38) 

where s is the unit vector in the direction of B. The Lorentz force then becomes 

Fl = --^VS2 + -1ss-V52 + :^s-Vs. (1.39) 
ovr ovr 4vr 

Combining the first two terms formally, we can write this as 

Fl = -^V^B^ + • Vs, (1.40) 
ovr 4vr 
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where V_l is the projection of the gradient operator on a plane perpendicular to B. The 
first term, perpendicular to the field lines, now describes the action of magnetic pressure 
more accurately. The second term, also perpendicular to B contains the effects of field 
line curvature. Its magnitude is 

B^ 

I— s-Vs| = -— , (1.41) 

where 

i?^ = l/|s-Vs| (1.42) 
is the radius of curvature of the path s. [problem 3.8: magnetic forces in a 1/r^ field] 



As an example of magnetic curvature forces, consider an axisymmetric azimuthally 
directed field, 'B = B (j), in cylindrical coordinates (ro, 0, z). The strength B is then a 



function of zu and z only (see problem 3.9 1). The unit vector in the azimuthal direction 



(f) has the property cf) ■'^(f) = —tSj/vc!, so that 

-^(B- V)B = --^— ^. (1.43) 
At: A-K w 

The radius of curvature of the field line is thus the cylindrical radius vd. The curvature 
force is directed inward, toward the center of curvature. For an azimuthal field like this, 
it is often also referred to as the hoop stress (like the stress in hoops keeping a barrel 



together), [problem 3.9 magnetic forces in a 1/m field.] 



1.3.2. Magnetic stress tensor 

The most useful alternative form of the Lorentz force is in terms of the magnetic stress 



tensor. Writing the vector operators in terms of the permutation symbol e (sect. 4.1), 
one has 

[(V X B)xB]i = ^ijk^jim^^Bk 

dxi 

= ^^{BiBk-^B^ik), (1.44) 

where the summing convention over repeated indices is used and in the last line V • B = 
has been used. Define the magnetic stress tensor M by its components: 



Mi,- = —B'^6ii - ^BiBi. (1.45) 
^ Svr ^ 47r ^ ^ ^ 



Eq. (1.44) then shows that the force per unit volume exerted by the magnetic field is 

3r: 

(V X B)xB = -V • M. (1.46) 



minus the divergence of this tensor: 

1 



An 

The magnetic stress tensor thus plays a role like that of the fluid pressure in ordinary 
fluid mechanics, except that it is a rank -2 tensor instead of a scalar. This is much like 
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Figure 1.6.: Magnetic force vectors exerted by a uniform magnetic field 



the stress tensor in the theory of elasticity. If F is a volume bounded by a closed surface 
S, (1.46) yields by the divergence theorem 



j |^(V X B)xB dy = ^ -n-M dS, 



(1.47) 



where n is the outward normal to the surface S. This shows how the net Lorentz force 
acting on a volume V of fluid can be written as an integral of a magnetic stress vector, 



the integrand of the right hand side in (1.46), acting on its surface. If, instead, we are 
interested in the forces Fg exerted by the volume on its surroundings, a minus sign is to 
be added, 



Fs = n-M = —B^n - —BBn, (1.48) 
* Svr 47r ' ^ ^ 

where = B • n is the component of B along the outward normal n to the surface of 
the volume. 

The stress in a magnetic field differs strongly from that in a fluid under pressure. 
Unlike normal fluid pressure, magnetic stress does not act perpendicular to a surface 
but is a vector at some angle to the normal of the surface, as in a sheared elastic 



medium or viscous fluid (problem 3.10). 



It is often useful to visualize magnetic forces in terms of the surface stress vector ( 1 .48 ) , 



evaluated on surfaces of suitable orientation. The sign of this surface force depends on 
the direction of the normal taken on the surface. The ambiguity caused by this is removed 
by deciding on the physically relevant volume to which the surface element would belong 



(for an example see 2.11). 



1.3.3. Properties of the magnetic stress. Pressure and tension 

To get some idea of the behavior of magnetic stresses, take the simple case of a uniform 
magnetic field, in the z-direction say, and evaluate the forces exerted by a volume of this 
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magnetic field on its surroundings, see Fig. 1.6 The force F5 whicli tlie box exerts on 



a surface parallel to B, for example the surface perpendicular to the x-axis on the right 



side of the box (1.48 with n = x), is Fright = x • M. The components are 



-Pright.x — ^-^^ ~ in'^^^^ ~ Svr'^^' -^"ght,2 — -bright,?/ — 0. (1.49) 

Only the magnetic pressure term contributes on this surface. The magnetic field exerts 
a force in the positive j;-direction, away from the volume. The stress exerted by the 
magnetic field at the top surface of the box has the components 



^top,2 = ^-^^ ~ ^^z^z — Ptop,x = Ftop,y = 0, (1.50) 



i.e. the stress vector is also perpendicular to the top surface. It is of equal magnitude to 
that of the magnetic pressure exerted at the vertical surfaces, but of opposite sign. 

On its own, the magnetic pressure would make the volume of magnetic field expand in 
the perpendicular directions x and y. But in the direction along a magnetic field line the 
volume would contract. Along the field lines the magnetic stress thus acts like a negative 
pressure, as in a stretched elastic wire. As in the theory of elasticity, this negative stress 
is referred to as the tension along the magnetic field lines. 

A magnetic field in a conducting fluid thus acts somewhat like a deformable, elastic 
medium. Unlike a usual elastic medium, however, it is always under compression in two 
directions (perpendicular to the field) and under tension in the third (along the field 
lines), irrespective of the deformation. Also unlike elastic wires, magnetic field lines 
have no 'ends' and cannot be broken. As a consequence, the contraction of the box in 



Fig. |1.6| under magnetic stress does not happen in practice, since the tension at its top 
and bottom surfaces is balanced by the tension in the magnetic lines continuing above 
and below the box. The effects of tension in a magnetic field manifest themselves more 
indirectly, through the curvature of field lines (see eq. 1.40| ff). For an example, see 
sect.ETH 

Summarizing, the stress tensor plays a role analogous to a scalar pressure like the gas 



pressure, but unlike gas pressure is extremely anisotropic. The first term in (1.45) acts 



in the same way as a hydrodynamic pressure, but it is never alone. Approximating the 



effect of a magnetic field by a scalar pressure term is rarely useful (see also sect. 2.12). 



1.3.4. Boundaries between regions of different field strength 

Let be a surface in the field configuration, everywhere parallel to the magnetic field 
lines but otherwise of arbitrary shape (a so-called 'magnetic surface'). Introduce a 'total' 
stress tensor Mt by combining gas pressure with the magnetic pressure term: 

Mt,ij = ip + B^/^Ti)5,j - BiBj/An. (1.51) 

The equation of motion is then 

dv 

= -divMt + pg. (1.52) 

Integrate this equation over the volume of a 'pill box' V of infinitesimal thickness e 
which includes a unit of surface area of S (Fig. 1.7). The left hand side and the gravity 
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Figure 1.7.: Boundary between two regions of different field strength, with the pill box (red) 
used for applying Gauss's theorem. 



term do not contribute to the integral since the volume of the box vanishes for e — t- 0, 
hence 

divMtdy = 0. (1.53) 

pill 



The stress ( 1.51 ) is perpendicular to S on the surfaces of the box that are parallel to B. 
Applying Gauss's theorem to the box then yields, in the limit e — ?■ 0, 



Pi + 



Svr 



P2 + 



R2 

111 

Svr 



(1.54) 



Equilibrium at the boundary between the two regions 1 and 2 is thus governed by the 
total pressure p+B'^/Sn only. The curvature force does not enter in the balance between 
regions of different field strength (flag this as being contrary to your intuition). 



The configuration need not be static for (1.54) to apply, but singular accelerations at 



the surface S must be excluded. The analysis can be extended to include the possibility 
of a sudden change of velocity across S (so the contribution from the left hand side of 
1.52 does not vanish any more), and S can then also be taken at an arbitrary angle to 



the magnetic field. This leads to the theory of MHD shock waves. 



1.3.5. Magnetic buoyancy 

Consider a magnetic flux bundle embedded in a nonmagnetic plasma, in pressure equi- 



librium with it according to (1.54). With Be = 0, the pressure inside the bundle is thus 
lower than outside. Assume that the plasma has an equation of state p = TZpT, where 
TZ (erg g^^ K^^) is the gas constant and T the temperature, and assume that thermal 
diffusion has equalized the temperature inside the bundle to the external temperature 
T. The reduced pressure then means that the density pi is reduced: 



or 



SP= Pi- Pe 



6p/pe = -1/(3 



B^ 



1 /^Mx; n2 

2^ Ci ^ ' 



(1.55) 



(1.56) 



where Cj the isothermal sound speed q = {p/pY^"^ and v^e = B / {inpcW^'^ is a notional 
Alfven speed, based on the internal field strength and the external density. In the 
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presence of an acceleration of gravity g, the reduced density causes a buoyancy force Fb 
against the direction of gravity. Per unit volume: 



Fb = g5p. (1.57) 
This causes a tendency for magnetic fields in stars to drift outwards (see problems 3.13 



3.14 for the expected speed of this process). 



1.4. Strong fields and weak fields, plasma-/3 

In the examples so far we have looked separately at the effect of a given velocity field on 
a magnetic field, and at the magnetic forces on their own. Before including both, it us 
useful to classify physical parameter regimes by considering the relative importance of 
the terms in the equation of motion. Ignoring viscosity and external forces like gravity. 



the equation of motion (1.18) is 



dv 1 
p— = -Vp+ — (V X B)xB . (1.58) 
dt Att 

A systematic procedure for estimating the relative magnitude of the terms is to decide 
on a length scale / and a time scale r that are characteristic for the problem at hand, 
as well as characteristic values vq for velocity and Bq for field strength. Since the sound 
speed is generally a relevant quantity, we assume a compressible medium, for simplicity 
with an equation of state as before, p = TZpT, where TZ is the gas constant and T the 
temperature, which we take to be constant here ('isothermal' equation of state). Define 
dimensionless variables and denote them with a tilde ~: 

t = Tt, V = V/l, V = vqv, B = BqB. (1.59) 

With the isothermal sound speed ci, 

cf = p/p = TZT, (1.60) 
the equation of motion becomes, after multiplication by l/p, 

i;o--iv = -cfVlnp + ?ji(VxB)xB, (1.61) 
T dt 



where is the Alfven speed: 
A characteristic velocity for things happening at the time scale r over the length of 



interest I is vq = 1/t. Dividing by (1.61) by cf then yields a dimensionless form of the 
equation of motion: 

A^^iLv = -Vln/)+ |(VxB)xB, (1.63) 

dt p 

where M is the Mach number of the flow, 

M = vq/cu (1.64) 
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and (3 is the so-calle(^ plasma- P: 



the ratio of gas pressure to magnetic pressure. Since we have assumed that /, r, vq, Bq are 



representative values, the tilded quantities in (1.63) are all of order unity. The relative 
importance of the 3 terms in the equation of motion is thus determined by the values of 
the two dimensionless parameters M and (3. 



Assume first the case of a highly subsonic flow, 7W ^ 1, so the left hand side of (1.63) 
can be ignored. The character of the problem is then decided by the value of f3. The 
physics of high-(5 and low-(3 environments is very different. 

If /3 S> 1, i.e. if the gas pressure is much larger than the magnetic energy density, the 
second term on the right is small, hence the first must also be small, V In /? ^ 1. That is 
to say, the changes in density produced by the magnetic forces are small. In the absence 
of other forces causing density gradients, a constant-density approximation is therefore 
often useful in high-/3 environments. 

If /3 <C 1, on the other hand, the second term is large, but since logarithms are not 
very large, it cannot be balanced by the first term. We conclude that in a low-/?, low-A^ 
plasma the magnetic forces must be small: we must have 

(VxB)xB ~ ©(Z?) < 1. (1.66) 

In the limit /3 — )• this can be satisfied in two ways. Either the current vanishes, 
V X B = 0, or it is parallel to B. In the first case, the field is called a potential field. It 
has a scalar potential (j)^^ such that B = — V</)m, and with V - B = 0, it is governed by 
the Laplace equation 

VVm = 0. (1.67) 

Applications could be the magnetic field in the atmospheres of stars, for example. 
For reasons that are less immediately evident at this point, the currents in such an 
atmosphere are small enough for a potential field to be a useful first approximation, 
depending on the physical question of interest. The more general second case, 

(VxB)xB = (1.68) 



describes force-free fields (sect. 1.5.1). Like potential fields, they require a nearby 'an- 
choring' surface (sect. 1.5.3). They are also restricted to environments such as the 
tenuous atmospheres of stars like the Sun, magnetic A stars, pulsars, and accretion 
disks. 

If the Mach number M is not negligible and (3 is large, the second term on the right 



of (1.63) can be ignored. The balance is then between pressure forces and accelerations 
of the flow, i.e. we have ordinary hydrodynamics with the magnetic field playing only 
a passive role. We can nevertheless be interested to see how a magnetic field develops 
under the infiuence of such a flow. This is the study of the kinematics of a magnetic 
fleld, or equivalently: the study of the induction equation for different kinds of specified 
flows. 



^ The plasma-/3 has become standard usage, thanks to the early plasma physics literature. Attempts 
to introduce its inverse as a more logical measure of the influence of a magnetic field have not been 
very successful. 
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Finally, if (3 is small and the field is not a force-free or potential field, a balance 
is possible only if M is of order ^ 1. That is, the flows must be supersonic, 
with velocities vq ~ wa- The magnetic fields in star- forming clouds are believed to be 
approximately in this regime. 

The intermediate case /3 ~ 1 is sometimes called 'equipartition', the 'equi' referring 
to the approximate equality of magnetic and thermal energy densities in this case. The 
term equipartition is not unique, however; it is also used for cases where the magnetic 
energy density is comparable with the kinetic energy density of the flow, i.e. when 

or equivalently 

V « VA- (1.70) 



1.5. Force-free fields and potential fields 
1.5.1. Force-free fields 

In a force-free field, VxB is parallel to B. Hence there is a scalar a such that 

VxB = aB. (1.71) 

Taking the divergence we find: 

B-Va = 0, (1.72) 

that is, a is constant along field lines. Force-free fields are 'twisted': the field lines in 
the neighborhood of a field line with a given value of a 'wrap around' it, at a pitch 
proportional to a. 

The case of a constant a everywhere would be mathematically interesting, since it 
leads to a tractable equation. This special case (sometimes called a 'linear' force-free 
field) is of little use, however. Where force-free fields develop in nature, the opposite 
tends to be the case: the scalar a varies strongly between neighboring magnetic surfaces 



(cf. sects. 1.11, 2.4) 



With magnetic forces vanishing, infinitesimal fluid displacements ^ do not do work 
against the magnetic field, hence the energy in a force-free magnetic field is an extremum. 
The opposite is also the case: if the magnetic energy of a configuration is a (possibly 
local) minimum, the field must be force-free. 

To make this a bit more precise imagine a (closed, simply connected) volume V of 
perfectly conducting fluid in which displacements ^ take place, surrounded by a volume in 
which the fluid is kept at rest, so that ^ = there. We start with an arbitrary magnetic 
fleld configuration in V, and let it relax under this constraint until an equilibrium is 
reached. This relaxing can be done for example by adding a source of viscosity so fluid 
motions are damped out. In doing so the field lines move to different locations, except at 
their ends, where they are kept in place by the external volume. To be shown is that the 
minimum energy configuration reached in this way is a force- free fleld, (VxB)xB = 0. 

Small displacements inside the volume change the magnetic energy £"111 in it by an 
amount 5E^, and the magnetic field by an amount 5B: 



Ait 6E„ 



B^dV 



B • 5B dV. (1.73) 
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The velocity of the fluid is related to the displacements ^ by 

V = di/dt, (1.74) 
so that for small displacements the induction equation is equivalent to 

5B = Vx(^xB). (1.75) 



(problem 3.5). Using one of the vector identities and the divergence theorem, (1.73) 
yields 

47r 5E^ = Jb-Vx{$xB)= (1.76) 

= X B) X B] • dS + y"(V X B) • (^ X B) dV. (1.77) 
The surface term vanishes since the external volume is being kept at rest, so that 

47r6E^ = y"(Vx B) • (^ X B) dy = -y"^-[(Vx B)xB] dF. (1.78) 

If the energy is at a minimum, 6E^ must vanishes to first order for arbitrary displace- 
ments This is possible only if the factor in square brackets vanishes everywhere, 

(VxB)xB = 0, (1.79) 

showing that a minimum energy state is indeed force-free. To show that a force-free field 
in a given volume is a minimum, rather than a maximum, requires examination of the 
second order variation of 5E^ with ^ (cf. Roberts 1967, Kulsrud 2005)[^ 



1.5.2. Potential fields 

The energy of a force- free minimum energy state can be reduced further only by relaxing 
the constraint of perfect conductivity. Assume that the magnetic field in the external 
volume is again kept fixed, for example in a perfectly conducting medium. Inside V the 
lines can slip with respect to the fiuid, by allowing magnetic diffusion to take place. The 
only constraint on the magnetic field inside V is now that it is divergence-free. We can 
take this into account by writing the changes in terms of a vector potential, 

5B = Vx 5A. (1.80) 

The variation in energy, 5E^ then is 

inSE^ = [ B-Vx6AdV (1.81) 
Jv 

= [ V ■ {5AxB) dV + [ <5A-VxBdy (1.82) 
Jv Jv 

= (f SA X B dS + I 5A-VxBdV. (1.83) 
Js Jv 

* In addition there is a question of uniqueness of the minimum energy state. This question is not quite 
trivial (cf. Moffatt 1985). 
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In order to translate the boundary condition into a condition on 5 A, a gauge is needed. 



Eq. (1.75) shows that in the perfectly conducting external volume, we can take 6 A to 
be 

6A = ^xB. (1.84) 

In this gauge the condition that field lines are kept in place on the boundary, ^ x B = 0, 
thus requires setting 6 A = on the boundary. The surface term then vanishes. Since 
we have no further constraints, 5 A is otherwise arbitrary inside V, so that the condition 
of vanishing energy variation is satisfied only when 

VxB = (1.85) 

inside V. That is, the field has a potential (pj^, B = — V(/)m, governed by the Laplace 
equation, 

VVm = 0. (1.86) 

The boundary condition on (p is found from V • B = 0. The component of B normal 
to the surface n, n • B = d(p/dn, is continuous across it. From potential theory it then 
follows that is unique. There is only one potential for such Neumann- type boundary 
conditions, up to an arbitrary constant. This constant only affects the potential, not the 
magnetic field itself. The minimum energy state of a magnetic field with field lines kept 
fixed at the boundary is thus a uniquely defined potential field. 

Potential fields as energy source? 

A consequence of the above is that the magnetic energy density B'^ /Sn of a potential field 
is not directly available for other purposes. There is no local process that can extract 
energy from the lowest energy state, the potential field. Its energy can only be changed 
or exploited by changes in the boundary conditions. If the fluid is perfectly conducting, 
the same applies to a force-free field. 

1.5.3. The role of the boundaries in a force-free field 

The external volume that keeps the field fixed in the above is more than a mathematical 
device. Its presence has a physical significance: the magnetic field exerts a stress on the 



boundary surface, as discussed in section p^. 3. 2[ The external medium has to be able to 
take up this stress. 

Irrespective of its internal construction, a magnetic field represents a positive energy 
density, making it expansive by nature. The internal forces (the divergence of the stress 
tensor) vanish in a force- free field, but the stress tensor itself does not; it vanishes only 
where the field itself is zero. At some point there must be something else to take up 
internal stress, to keep a field together. In the laboratory this is the set of external 
current carrying coils. In astrophysics, the magnetic stress can be supported by, for 
example a stellar interior, a gravitating cloud or by the surface of an accretion disk. The 
intrinsically expansive nature of magnetic fields can be formalized a bit with an equation 
for the global balance between various forms of energy, the tensor virial equation (e.g. 
Kulsrud Chapter 4.6). 

Whereas a potential field is determined uniquely by the instantaneous value of the 
(normal component of) the field on its boundary, the shape of a force-free field con- 
figuration also depends on the history of things happening on its boundary. Rotating 
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displacements on the boundary wrap the field lines inside the volume around each other. 
The values of a which measure this wrapping reflect the entire history of fluid displace- 
ments on the boundary surface. 

Since the value of a is constant along a field line, it is also the same at its points of 
entry and exit on the boundary. The consequence is that, unlike in the case of a potential 
field, the construction of a force-free field is not possible in terms of a boundary-value 
problem. A given force-free field has a unique distribution of a and Bn on its surface. But 
there is no useful inverse of this fact, because the correspondence of the points of entry 
and exit of the field lines is not known until the force- free field has been constructed. 
Force-free fields must be understood in terms of the history of the fluid displacements 



at their boundary. See section 2.4.1 for an application. 



1.5.4. The vanishing force-free field theorem 

A consequence of the expansive nature of magnetic fields is the following Theorem: A 
force-free field which does not exert stress at its boundaries vanishes everywhere inside. 
To show this examine the following integral over a closed volume V with surface S: 



^ J x,[(VxB)xB], dF = J 



Xj ^ dV 



d 

dxk 



(xiM. 



dV 



dxk 

XiMjkUk dS — I Mij dV, 



(1.87) 



where are again the components of the outward normal to the surface S. For a 
force-free field, Oij = 0. Taking the trace of (1.87) we find 



XiMikUk dS= I MiidV = ^ I dV. 



(1- 



If the stress Mj^ vanishes everywhere on the surface 5 (left hand side), it follows that 
-6 = everywhere inside V . 

Taking the surface of V to infinity, the theorem also implies that it is not possible to 



construct a field that is force-free everywhere in an unbounded volume. (Problem 3.17 
models for magnetic A stars). 

Since a force-free field has a higher energy density than a potential field, for given 
boundary conditions, twisting a magnetic field does not help to 'keep it together', con- 
trary to possible expectation. One might think that the 'hoop stresses' caused by twisting 
might help, much like an elastic band can be used to keep a bundle of sticks together. 
Instead, the increased magnetic pressure due to the twisting more than compensates the 
hoop stress. This is illustrated further in the next section. 

Summarizing: like a fiuid under pressure, a magnetic field has an internal stress M. 
The divergence of this stress may vanish, as it does in a force-free field, but there still 
has to be a boundary capable of taking up the stress exerted on it by the field. If the 
stress cannot be supported by a boundary, it has to be supported by something else 
inside the volume, and the field cannot be force-free. 



Problem 3.11 stress on the surface of a uniformly magnetized sphere. 
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Figure 1.8.: Bundle of field lines ('flux tube') embedded in pressure equilibrium in a field-free 
plasma. Twisting it (right) by applying a torque (arrows) at top and bottom causes 
its width to expand. The net current along the tube vanishes. 



1.6. Twisted magnetic fields 



1.6.1. Twisted fluxtubes 



Fig. 1.8 shows a bundle of field lines embedded in a nonmagnetic plasma, with field 
strength Bq = Bqz along the z-axis. If pi and pe are the internal and external gas 
pressures, the condition for the tube to be in pressure equilibrium is pi + Bq/8tt = p^ 



(see 1.3.4). We twist a section of the tube by rotating it in opposite directions at top 



and bottom. Imagine drawing a closed contour around the tube in the field-free plasma. 



With (1.13) and applying Stokes' theorem to this contour: 



dS 



B • dl = 0, 



(1.89) 



where / is the path along the contour, and S a surface bounded by the contour. The 
net current along the tube vanishes no matter how the tube is twisted (see also problem 



3.16) 



The twisting has produced an azimuthal component Bfj, = Sgtanx, where x is the 
pitch angle of the twisted field. The magnetic pressure at the boundary of the tube 
has increased: = Bq + so the tube is not in pressure equilibrium anymore (see 



problem 3.15). The additional pressure exerted by B^ causes the tube to expand. 



One's expectation might have been that the tube radius would contract due to the 
hoop stress of the added component B^. We see that the opposite is the case. A 



view commonly encountered is that the twist in Fig. 1.8 corresponds to a current flowing 



along the axis, and that such a current must lead to contraction because parallel currents 
attract each other. The hoop stress can in fact cause some of the field configuration in 
the tube to contract (depending on how the twisting has been applied) , but its boundary 



always expands unless the external pressure is also increased (see problem 3.16 to resolve 
the apparent contradiction). This shows how thinking in terms of currents as in a 
laboratory setup leads astray. For more on this and the popular mistakes made in this 
context see Ch. 9 in Parker (1979). For an illustration of the above with the example of 



MHD jets produced by rotating objects, see section 2.13 
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1.6.2. Magnetic helicity 



If A is a vector potential of B, the magnetic helicity H of a field configuration is defined 
as the volume integral 



H 



A-BdV. 



(1.90) 



Its value depends on the arbitrary gauge used for A. It becomes a more useful quantity 
when the condition 

B-n = (1.91) 

holds on the surface of the (simply connected) volume V. There are then no field lines 
sticking through its surface, the field is completely 'contained' within V. In this case H 
(dimensions G^cm'*), has a definite value independent of the gauge. [For proofs of this 
and related facts, see section 3.5 in Mestel (2012)]. It is a global measure of the degree 
of twisting of the field configuration. Because of the gauge dependence, H does not have 
locally definable values: there is no 'helicity density'. 

This is not just a computational inconvenience. The physically significant reason is 
that H is a topological quantity: it depends not only on a local degree of twisting, but 
also on the global 'linking' of field lines within the configuration. This is illustrated in 



Fig. 1.9 showing two loops of magnetic field, with magnetic fluxes <I>i and $2- If they 
are not linked (panel a), each of the loops can be fit inside its own volume satisfying 



( 1.91 ), hence the helicity is a sum over the two loops. If they are untwisted, in the sense 



that all field lines close on going once around the loop without making turns around the 
cross section of the tube, the helicity of each vanishes. If the loops are linked, however, 
(panel b) it can be shown that H = 2<I>i<I>2 (problem 3.21). For an intuitive picture 
how this comes about in terms of the definition ( |1.90 ), inspect the location of loop 2 
relative to the vector potential of loop 1 . In case (b) , the field of loop 2 runs in the same 
direction as the field lines of Ai, in case (a), much of it runs in the direction opposite 
to Ai. 

The importance of H lies in its approximate conservation property. In ideal MHD, H 
is an exactly conserved quantity of a magnetic configuration. Within a volume where 



(1.91) holds, H does not change as long as perfect flux freezing holds everywhere. In 



practice, this is not a very realistic requirement, since flux freezing is rarely exact over 
the entire volume. Reconnection (sects. 1.11| |2.4[ ) is bound to take place at some point 





Figure 1.9.: Helicity of configurations of two untwisted magnetic loops. Green: a field line of 
the vector potential Ai of the left loop. In the left panel the loops are not linked 
and the helicity of the configuration vanishes. When the loops are linked (b), the 
helicity does not vanish. 
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in time at some place in the volume, causing the topology of the field to change, or 



it may not be possible all the time to find a volume where (1.91) holds. It turns out, 
however, that H is often conserved in a more approximate sense. In the laboratory, a 
helical magnetic field configuration that is not in equilibrium, or unstable, first evolves 
on a fast time scale of a few Alfven crossing times L/va- During this phase H is 
approximately conserved. Following this, the magnetic helicity evolves more slowly by 
reconnection processes. 

Helicity does not behave like magnetic energy. In the presence of reconnection, H 
can decrease during a relaxation process in which magnetic energy is released, but it can 
also increase. The evolution of magnetic fields on the surface of magnetically active stars 
(like the Sun) is sometimes described in terms of 'helicity ejection'. In view of the global, 
topological nature of magnetic helicity, this usage of the term helicity is misleading. In 
a magnetic eruption process from the surface of a star magnetic helicity can decrease 



even when the thing being ejected is not helical at all. For more on this see section 2.14[ 



There are other ways of characterizing twist, for example the quantity h = (VxB) 



B/i? with dimension 1 /length, called current helicity. In a force- free field (sect. 1.5.1) 
its value is equal to a. In contrast to H which is a property of the configuration as 
a whole, current helicity is a locally defined quantity. It does not have a conservation 
property, however. 

Ultimately, these facts are all a consequence of the non-local (solenoidal vector-) nature 
of the magnetic field itself. The usefulness of analogies with the conservation properties 
of other fiuid quantities is intrinsically limited. 

1.7. Stream function 

Though in practice all magnetic fields in astrophysics are 3-dimensional in one or the 
other essential way, 2-dimensional models have played a major role in the development 
of astrophysical MHD, and their properties and nomenclature have become standard 
fare. Historically important applications are models for steady jets and stellar winds. 

If (ro, (p, z) are cylindrical coordinates as before, an axisymmetric field is independent 
of the azimuthal coordinate: 

dB/d4> = 0. (1.92) 

(In numerical simulations axisymmetric models are sometimes called '2.5-dimensionar). 
Such a field can be decomposed into its poloidal and toroidal components 

B = Bp(ro,z) + Bt(w,z), (1.93) 

where the poloidal field contains the components in a meridional plane (p = est.: 

Bp = (5^,0,5,), (1.94) 

and 

Bt = B^4, (1.95) 

is the azimuthal field component (the names toroidal and azimuthal are used inter- 
changeably in this context). Define the stream functioi^ a scalar tp, by 

ip{w,z)= I wBzdw. (1.96) 
Jo 



^ A stream function can be defined for any axisymmetric solenoidal vector. In fiuid flows, it is called 
the Stokes stream function. 
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Using divB = 0, the poloidal field can be written as 



= —dip/dw, = -—d^jj/dz. (1.97) 

TU W 

From this it follows that -0 is constant along field lines: B • V0 = 0. The value of ■0 can 
be used to label a field line (more accurately: an axisymmetric magnetic surface). It 
equals (modulo a factor 27r) the magnetic flux contained within a circle of radius w from 
the axis. It can also be written in terms of a suitable axisymmetric vector potential A 
of B: 

= wA^, (1.98) 



but is a more useful quantity for 2D configuration^^ Stream functions can also be 



defined in planar 2D symmetry (problem 3.22). 

1.8. Waves 

A compressible magnetic fluid supports three types of waves. Only one of these is 
comparable to the sound wave familiar from ordinary hydrodynamics. It takes some time 
to develop a feel for the other two. The basic properties of the waves modes of a uniform 
magnetic fluid also turn up in other MHD problems, such as the various instabilities. 
Familiarity with these properties is important for the physical understanding of time- 
dependent problems in general. 

The most important properties of the waves are found by considering first the simplest 
case, a homogeneous magnetic field B in a uniform fluid initially at rest (v = 0). The 
magnetic fleld vector defines a preferred direction, but the two directions perpendicular 
to B are equivalent, so the wave problem is effectively two-dimensional. In Cartesian 
coordinates {x,y,z), take the initial magnetic field (also called 'background' magnetic 
field) B along the z-axis, 

B = Bz (1.99) 

where i? is a constant. Then the x and y coordinates are equivalent, and one of them, 
say y, can be ignored by restricting attention to perturbations 6q that are independent 
of y: 

dy5q = 0. (1.100) 

Write the magnetic field as B-|-(5B, the density as p + 5p, the pressure as p + 5p, where 
5B, 6p and 5p are small perturbations. Expanding to first order in the small quantities 
6q (of which v is one), the linearized equations of motion and induction are: 

dv 1 
p— = -V6p+ —{V X 5B)xB. (1.101) 

at 47r 



dt 

The continuity equation becomes 



V X (v X B). (1.102) 



^ + pV-v = (1.103) 



Plotting the contour lines i/i = est., for example, is an elegant way to visualize field lines in 2-D. 
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(since the background density p is constant). To connect 5p and 5p, assume that the 
changes are adiabatic: 

h = h = ^2^/5, (1.104) 



where the derivative is taken at constant entropy, and Cs is the adiabatic sound speed. 
For an ideal gas with ratio of specific heats 7 = Cp/c„, 

cl = jp/p. (1.105) 

The components of the equations of motion and induction are 

dvx d6p B d6Bx d5Bz. dvz d6p 

86 B„ „dvn, 

—l = B-^. 1.109 
ot az 

The y-components only involve 5By and Vy. As a result, eqs. (1.103 - 1.109) have 



solutions in which Vx = Vz = 6Bx = 5Bz = 5p = 6p = 0, with dBy and Vy determined by 
(1.107) and (1.109). These can be combined into the wave equation 

^-^A^)mj^vy) = 0, (1.110) 
where 

and the amplitudes are related by 

6By/B = ± Vy/VA- (1.112) 

These solutions are called Alfven waves (or 'intermediate wave' by some authors). Since 
they involve only the y and z coordinates, they 'propagate in the y — z plane only'. 

The second set of solutions involves Vx,Vz, 6Bx, SBz, 6p, 6p, while Vy = 5By = 0. Since 
the undisturbed medium is homogeneous and time- independent, the perturbations can 
be decomposed into plane waves, which we represent in the usual way in terms of a 
complex amplitude. Any physical quantity q thus varies in space and time as 

(7 = go exp[i(u;t — k • x)], (1.113) 

where u is the circular frequency and qq is a (complex) constant. The direction of the 
wave vector k (taken to be real) is called the 'direction of propagation' of the wave. With 
this representation time and spatial derivatives are replaced by iuj and — zk respectively. 
The continuity equation yields 

6p- pk-v/u = 0. (1-114) 
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With the adiabatic relation (1.104), this yields the pressure: 

6p = c^p k • v/u!. 



(1.115) 



Substitution in (1.106, 1.108) yields a homogenous system of linear algebraic equations: 

upvx = kxp{kxVx + kzVz)cl/uj - -—{kz5Bx - k^dBz) (1.116) 



47r 

(^PVz = kzPikxVx + kzVz)cl/uj 
uj5Bx = -kzBvx 
ojSBz = kxBvx- 



(1.117) 
(1.118) 
(1.119) 



Writing this as a matrix equation Aijqj = 0, where q = {vx,Vz,5Bx,SBz), the condi- 
tion that nontrivial solutions exist is det A = 0. This yields the dispersion relation for 
the compressive modes (the relation between frequency uj and wavenumber k): 



UJ 



oj^k\ct + vi) + klk^civi = 0. 



(1.120) 



Let 9 be the angle between k and B: 

kx = ksinO, 



kz = k cos 9, 



(1.121) 



and set u = u/k, the phase velocity of the mode. Then (1.120) can be written as 

n*^ - (c^ + f a) + c^vl cos^ 9 = 0. ( 



1.122) 



The four roots of this equation are real and describe the magnetoacoustic waves. There 
are two wave modes, the slow and the fast magnetoacoustic modes, or 'slow mode' and 
'fast mode' for short, each in two opposite directions of propagation (sign(a;//c) = ±1). 
They are also collectively called the magnetosonic modes. 



1.8.1. Properties of the Alfven wave 



With (1.113), (1.110) becomes: 



7 2 2 

kz^A, 



(1.123) 



the dispersion relation of the Alfven wave. Since this only involves kz one says that 
Alfven waves 'propagate along the magnetic field': their frequency depends only on the 
wavenumber component along B. This does not mean that the wave travels along a 
single field line. In the plane geometry used here, an entire plane x =cst. moves back 
and forth in the y-direction (Fig. 1.10). Each plane x = est. moves independently of the 



others (see problem 3.23). The direction of propagation of the wave energy is given by 
the group velocity: 

^ = ^Az. (1.124) 

That is, the energy propagates along field lines, independ ent of the 'direction of propa- 
gation' k (but again: not along a single field line). From (1.112) it follows that 



Sir 2 



(1.125) 
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Figure 1.10.: Linear Alfven wave propagating along a plane magnetic surface (wave amplitude 
6B/B exaggerated). 



so there is 'equipartition' between magnetic and kinetic energy in an Alfven wave (as in 
any harmonic oscillator). Further properties of the Alfven wave are: 

- A linear Alfven wave is incompressive, i.e. 5p = 0, in contrast with the remaining modes, 

- it is transversal: the amplitudes (5B and v are perpendicular to the direction of prop- 
agation (as well as being perpendicular to B), 

- it is nondispersive: all frequencies propagate at the same speed. 

- in an incompressible fluid, it remains linear at arbitrarily high amplitude (e.g. Roberts 
1967). 

An example of a plane linear Alfven wave traveling on a single magnetic surface is 



shown in Fig. 1.10 The amplitude of the wave has been exaggerated in this sketch. The 
magnetic pressure due to the wave, {6B)'^/87r excites a compressive wave perpendicular 
to the plane of the wave ('mode coupling'). The wave will remain linear only when 
6B/B < 1. 

An Alfven wave can propagate along a magnetic surface of any shape, not necessarily 
plane. For it to propagate as derived above, however, conditions on this plane have to 
be uniform (B,/j, Cg constant). Though it cannot strictly speaking travel along a single 
field line, it can travel along a (narrow) magnetic surface surrounding a given field line: 
a 'torsional' Alfven wave. 



Torsional Alfven waves 

Instead of moving an infinite plane x =cst., we can also produce a more localized Alfven 



wave by 'rotating a bundle of field lines'. This is illustrated in Fig. 1.11 At time t = 0, a 
circular disk perpendicular to the (initially uniform) magnetic field is put into rotation 
at a constant angular velocity 0,. A wave front moves up at the Alfven speed along 
the rotating field bundle. This setup corresponds to an Alfven wave of zero frequency: 
apart from the wave front it is time-independent. Of course, an arbitrary superposition 
of wave frequencies is also possible. Since the wave is nondispersive, the wave front 



remains sharp. See problems 3.25 and 3.24 
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Figure 1.11.: Torsional Alfven wave propagating along a cylindrical magnetic surface. 
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Figure 1.12.: Propagation diagram for the MHD waves in a uniform medium, showing the phase 
speed (length of the vector u) as a function of the angle 9 of the wave vector with 
respect to the direction of the magnetic field. 



1.8.2. Properties of the magnetoacoustic waves 

In the second set of waves, the density and pressure perturbations do not vanish, so 



they share some properties with sound waves. The phase speed u (1.122) depends on 
the sound speed, the Alfven speed and the angle between the wave vector and the 
magnetic field. Introducing a dimensionless phase speed u, 



U = tt(CsVA)^^^, 



(1.122) becomes 



Cs 



(1.126) 



(1.127) 



This shows that the properties of the wave can be characterized by just two parameters: 
the angle of propa gation 9 an d the ratio of sound speed to Alfven speed, c^/va = 
(7/3/2)^/^ (cf. eqs. 1.65, 1.105). The phase speed does not depend on frequency: like 
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Figure 1.13.: Propagation diagrams of the fast mode (left) and slow mode (right) for t;A = Cs- 
Curves (u) show the phase speed as in Fig. |1.12[ Arrows: direction of the group 
speed Vg as a function of the direction (dotted) of the phase speed. 



the Alfven wave, the waves are nondispersive. Since they are anisotropic, however, the 
phase speed is not the same as the group speed. The solutions of eq. 1.122) are 

1 



u 



0' 



-(c^ + fi)[l±(l-4cos^ 9/b 



2^1/21 



where 



Cs VA 

— + — > 2. 

VA Cs 



(1.128) 



(1.129) 



Hence is positive, and the wave frequencies real as expected (for real wave numbers 
k). The +(— ) sign corresponds to the fast (slow) mode. The limiting forms for b ^ 1 
(either because va 3> Cg or Cs ^ va) are of interest. In these limits the fast mode speed 
is 

U{ ^ {cl + viy/^, (cs > i;A or > Cs) (1.130) 

i.e. Uf is the largest of va and Cs, and is independent of 9: the fast mode propagates 
isotropically in these limiting cases. For va ^ Cs or Cs ^ va the slow mode speed 
becomes: 



.2„,2 



C^V 



d + vi 



cos^ 9, (cs ^» VA or VA ^ Cg 



(1.131) 



i.e. Us is smaller than both va and Cg in these limiting cases, and its angular dependence 
cos^ 9 is the same as that of the Alfven wave. 

Another interesting limiting case is cos^ 9 <C 1, that is, for wave vectors nearly per- 
pendicular to B. Expressions (1.130) and (1.131) hold in this limit as well, but they now 
apply for arbitrary va, Cg. In this context, (1.130) is called the fast magnetosonic speed 
while the quantity 

Cc = (1.132) 



(c2+t;2)l/2 



is called the cusp speed. 
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Va/Cs= 0.5 




Va/Cs=2.0 




Figure 1.14.: Directions of fluid displacement (red) in the fast mode (top) and slow mode (bot- 
tom) at low (lef t) an d high (right) va/cs- Grey: angular dependence of the phase 
speed as in Fig. 



1.12 



The nomenclature used in describing wave propagation can be a bit confusing in 
the case of the Alfven and slow waves. Because of their anisotropy, the wave vector 
k is not the direction in which the energy of the wave flows. The case cos^ <C 1 is 
called 'propagation perpendicular to B', though in the Alfven and the slow mode the 
wave energy actually propagates along B in this case. Propagation diagrams, shown in 



Fig. 1.12 can help visualization. These polar diagrams show the absolute value of the 



phase speed, as a function of the angle of k with respect to B. 

The curves for the Alfven wave are circles. The slow mode is neither a circle nor an 
ellipse, though the difference becomes noticeable only when va and Cg are nearly equal. 
This is shown in Fig. 1.13, the propagation diagrams for va/cs = 1- It also shows the 
direction of the group speed, as a function of the angle 6 of the wave vector. At va/cs = 1 
the group speed shows the largest variations in direction, both in the fast and in the 
slow mode. At other values of va/cs the direction of the group speed of the slow mode is 
close to B for all directions of the wave vector, like in the Alfven wave, while the group 
speed of the fast mode is nearly longitudinal, as in a sound wave. 



The symmetry between Cg and va suggested by (1.127) and the propagation diagrams 



in Figs. 1.12 and 1.13 is a bit deceptive. The behavior of the waves is still rather different 
for /3 > 1 and /3 < 1. For /3 > 1 the fast mode behaves like a sound wave, modified a 



bit by the additional restoring force due to magnetic pressure (Fig. 1.14, top left). At 



low /3 (top right), it propagates at the Alfven speed, with fluid displacements nearly 
perpendicular to B as in the Alfven wave, but in contrast with the Alfven wave nearly 
isotropically. 
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The slow mode has the most complex properties of the MHD waves. At low /3 (high 
va/cs) it behaves like a 1-dimensional sound wave 'guided by the field'. The fluid ex- 
ecutes a 'sloshing' motion along field lines (Fig. 1.14, lower right). In this limit sound 
waves propagating along neighboring field lines behave independently of each other. The 
dominance of the magnetic field prevents the wave's pressure fluctuations from causing 
displacements perpendicular to the field (in other words the nonlinear coupling to other 
waves is weak). 

In the high-/3 (low v^/cs) limit, the propagation diagram of the slow wave approaches 
that of the Alfven wave. The fiuid displacements (Fig. [1.14 lower left) become perpen- 
dicular to the wave vector. The fiow is now in the x — z plane instead of the y-direction. 
At high /3, the slow mode can be loosely regarded (in the limit /3 — )• oo: exactly) as a 
second polarization state of the Alfven wave (polarization here in the 'wave' sense, not 
the polarization of a medium in an electric field as described in sect. 2.10). 



Summary of the magnetoacoustic properties 

The somewhat complex behavior of the modes is visualized most easily in terms of the 
asymptotic limits va ^ Cg and va ^ Cg'- 

- The energy of the fast mode propagates roughly isotropically; for va ^ Cg with the 
speed of an Alfven wave and displacements perpendicular to the field (Fig, l.lij top 



right), for va ^ Cg with the speed and fiuid displacements of a sound wave (top left). 
- The energy of the slow mode propagates roughly along field lines; for va ^ Cg at the 



sound speed with displacements along the field (Fig. 1.14, bottom right), for va ^ Cg at 
the Alfven speed and with displacements that vary from perpendicular to B to parallel, 
depending on the wave vector (bottom left). 



Waves in inhomogeneous fields 

The waves as presented in this section are of course rarely found in their pure forms. 
When the Alfven speed or the sound speed or both vary in space, the behavior of MHD 
waves is richer. The subject of such inhomogeneous MHD waves has enjoyed extensive 
applied-mathematical development. Important concepts in this context are 'resonant 
absorption' and 'linear mode conversion'. 



1.9. Poynting flux in MHD 

The Poynting flux S of electromagnetic energy: 

S = -^ExB, (1.133) 

is usually thought of in connection with electromagnetic waves in vacuum. It is in fact 
deflned quite generally, and has a nice MHD-speciflc interpretation. With the MHD 
expression for the electric held, E = — v x B/c, we have 

S = -'-B X (v X B). (1.134) 

S thus vanishes in flows parallel to B. Writing out the cross-products, and denoting by 
v± the components of v in the plane perpendicular to B: 

S = v_L — . (1.135) 
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An MHD flow does not have to be a wave of some kind for the notion of a Poynting 
flux to apply. It also applies in other time dependent flows, and even in steady flows (the 
MHD flow in steady magnetically accelerated jets from accretion disks or black holes for 



example). See problem 3.27 



Expression (1.135 ) suggests an interpretation in terms of magnetic energy being carried 
with the fluid flow. But the magnetic energy density is B^/Stt. To see where the missing 
factor 2 comes from, consider the following analogy from ordinary hydrodynamics. In 
adiabatic flow, i.e. in the absence of dissipative or energy loss processes so that entropy 
is constant in a frame comoving with the flow, the equation of motion can be written 
aO 

dv 

— = —\'w + (additional forces)/p, (1.136) 

where w is the heat function or enthalpy of the fluid, 

w = {e+p)/p, (1.137) 

e is the internal (thermal) energy per unit volume, and p the gas pressure. At constant 
entropy, w satifies = Wp/ p. Conservation of energy in a flow can be expressed in 
terms of the Bernoulli integral. In the absence of additional forces, it is given by: 

E=\v'^ + w. (1.138) 

In a steady adiabatic flow, E can be shown to be constant along a flow line (the path 
taken by a fluid element). 

The above shows that the thermodynamic measure of thermal energy in hydrodynamic 
flows is the enthalpy, not just the internal energy e. What is the physical meaning 
of the additional pressure term in the energy balance of a flow? In addition to the 
internal energy carried by the flow, the p work done at the source of the flow must be 
accounted for in an energy balance: this adds the additional p. The analogous measure 
of magnetic energy in MHD flows would be: 

Wm = em+Pm, (1.139) 

where Cm and pm are the magnetic energy density and magnetic pressure. Both are 



equal to S^/Svr, hence = B'^/Att. The Poynting flux in MHD, ( 1.135 ) can thus 
be interpreted as a flux of magnetic enthalpy in the plane perpendicular to B. Note, 
however, that it cannot be simply added to the hydrodynamic enthalpy (in a Bernoulli 
integral, for example), since the two flow in different directions. Like the Poynting flux 
itself, is not needed for the technical aspects of MHD calculations but is useful for 
the interpretation of the results. 

1.10. Magnetic difrusion 

Assume (without justification for the moment) that in the fluid frame the current is 
proportional to the electric field: 

j' = cTcE', (1.140) 



for the following see also Landau & Lifshitz §5 
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Figure 1.15.: Flow of fluid 'across field lines' when magnetic difi^usion is present. Though not 
carried with the flow, the fleld lines can still be labeled in a time-independent way 
if there is a surface (shaded) where diffusion is sufficiently slow for the field lines 
to be 'frozen-in' on the time scales of interest. 



that is, a linear 'Ohm's law' applies, with electrical conductivity cJc (not to be confused 
with a the charge density). In the non-relativistic limit (where B = B', j = j' = 
c/47r V X B): 

(1.141) 



— = E' = (E + - X B), 



and the induction equation (1.4) becomes 



aB 

~dt 



V X [v X B - ryV X B], 



where 



47rcJr ' 



(1.142) 



(1.143) 



with dimensions cm^/s is called the magnetic diffusivity. Eq. (1.142) has the character 
of a diffusion equation. If r] is constant in space: 



5iB = V X (v X B) + 7?V^B. (f] = est.) 



(1.144) 



The second term on the right has the same form as in diffusion of heat (by thermal 



conductivity), or momentum (by viscosity). Eq. (1.142) is a parabolic equation, which 



implies that disturbances can propagate at arbitrarily high speeds. It is therefore not 
valid relativistically; there is no simple relativistic generalization. 



In the presence of diffusion (77 7^ 0) Alfven's theorem (1.27) no longer applies. The 



flux of magnetic field lines through a loop comoving with the fluid (as in Fig. 1.2) is no 



longer constant in time. Individual field lines cannot be identified anymore with the 
fluid elements 'attached' to them, which was so useful for visualizing time dependent 
magnetic fields. 

Even when diffusion is important, however, it is often possible to label field lines as if 
they had an individual, time independent identity. It is sufficient that somewhere in the 
volume considered there is a region where where field lines can be associated with fluid 
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elements in a time-independent manner. This could for example be a solid conducting 
surface, or a fluid region of high conductivity. In such a region field lines can again be 
labeled as before. Thanks to divB = 0, this label can then be extended along the entire 
length of each field line. We can then speak of a speed at which field lines diffuse across 
a fluid (if the fluid is more or less at rest), or how fast a fluid flows across the field (if 
the field is more or less stationary). The latter case is illustrated in Fig. 
The importance of diffusion, i.e. the influence of the second term in 



1.15 



quantified with a dimensional analysis, in the same way as was done in sect. 1.4 for the 



1.142) can be 



equation of motion. If L and V are characteristic values for the length scale and velocity 



of the problem, this yields a nondimensional form of (1.144): 



§ = Vx(vxB) + -^V^B, 

Ot Km 



(1.145) 



where Rm is the magnetic Reynolds number: 



Rrn = LV/7J. 



(1.146) 



If Rm is large, the most common case in astrophysics (exercise 3.26), the diffusion 



term can be ignored and ideal MHD used except in subvolumes where small length 
scales develop, such as in current sheets. 



The form of the diffusion term in (1.142) is appropriate for processes that can be rep- 



resented by an 'Ohm's law'. Among other things, this assumes that the mean velocities 
of all components of the fluid (the neutrals and the two charge carriers) are the same to 
sufficient accuracy. This is not the case if the mean velocity of the electrons due to the 
current they carry, or of the neutral component relative to the charged components, is 



not negligible. In these cases additional terms appear in (1.142), corresponding to Hall 



drift and ambipolar diffusion, respectively (see 2.9). 



1.11. Current sheets 

When field lines of different directions get close together, the finite conductivity of the 
plasma has to be taken into account. Field lines are no longer tied to the fiuid, and 
magnetic energy is released. Normally the topology of the field lines also changes: they 
reconnect. 



A prototypical case is shown in Fig. 1.16 a current sheet separates magnetic field 



lines running in opposite directions. The sheet runs perpendicular to the plane of the 
drawing. The close-ups illustrate that significant differences can exist in its internal 
structure. In a pressure-supported sheet (top right) the field strength vanishes in the 
middle. Equilibrium of the structure is provided by a varying fiuid pressure p, such that 
p + B^/8tt = est across the sheet. The current is perpendicular to the magnetic field 
(and to the plane of the drawing). 

In a force-free current sheet on the other hand (botton right), fiuid pressure is constant 
(more generally: small, f3 <^ 1). The absolute value B of the field strength is constant 
across the sheet, and the current is parallel to the field lines instead of perpendicular. 
Instead of the field strength vanishing in the middle, the direction of the field lines rotates 
over 180°. In this case the space between the opposing field lines is filled with other field 
lines, instead of with fiuid. The mechanism of reconnection is correspondingly different 



in the two cases (sect. 2.4). 
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Figure 1.16.: The two types of current sheets at a 180° change in direction of the field (left). 

Top right: a two-dimensional sheet perpendicular to the plane of the drawing. 
B vanishes in the middle, the current is perpendicular to B and to the plane; 
current density indicated in green. Bottom right: a force- free sheet is intrinsically 
3-dimensional. The current is in the plane of the sheet, parallel to B, and B is 
constant across the layer. 
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2. Supplementary 



2.1. Alfven's theorem 

This somewhat intuitive derivation fohows Ferraro &: Plumpton. See also Kulsrud. Con- 



sider, as in section 1.2, a surface S{t) bounded by a loop of fluid elements C(t), which 
changes in shape and position with time. Denote the path length along C by s, the unit 
vector along C by s. The flux through the loop is 

m= [ B(r,t)-dS, (2.1) 

Js 

where S = Sn, with n a unit vector normal to S. To be shown is that $ is constant in 
time if the magnetic field satisfies the ideal MHD induction equation S^B = V x (v x B). 



Fig. 2.1 shows the loop at two instants in time separated by an infinitesimal interval 
6t. We are interested in the difference 5$ in the flux through the loop between these 
two instants : 

5$= [ B{r,t + 5t)-dS- [ B(r,t)-dS. (2.2) 

Jsit+5t) Js 

We write this as 

(5$ = /i + /2, (2.3) 



with 

h= [ B{r,t + 6t) ■ dS - [ B{r,t + 5t)-dS, (2.4) 

Js{t+st) Js{t) 

l2= [ B{r,t + St) ■ dS - [ B(r,t)-dS. (2.5) 

Js{t) Js{t) 

In this way we have divided the change in ^ into a part (/i) that is due to the change 
in position of the loop, and a contribution I2 due to the change of B in time at a fixed 
location. Hence 

h^StJ^-^iS. (2.6) 

We now apply the divergence theorem to the volume V defined by the surfaces S{t+St), 
S{t), and the side surface a traced out by the loop in the course of its displacement from 
t to t + 6t. The total flux through these surfaces is 



L 



B-dS- / B-dS+ / B-dS= / V-Bdr = 0. (2.7) 

S{t+St) Js{t) Jtj Jv 

The minus sign in the second term appears because we want the flux $ to have the same 
sign at the two instants t and t + 6t. With the sign of $ chosen to correspond to the 



outward normal on the upper surface in Fig. 2.1 and the outward normal being required 



for application of the theorem, a change of sign is thus needed for the lower surface, 
represented by the second term. 
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Figure 2.1.: Flux loops for heuristic derivation of Alfven's theorem. Red: flow lines 



To evaluate the flux through the side surface, we divide it into infinitesimal surface 
elements spanned by the vectors ds and vSt, where s is the arc length along the loop. 
The sign of ds is chosen such that the vector dS = ds x v6t has the direction of the 
outward normal on a. With these definitions we have 



/ B-dS = (f B-(ds X v)5t = 6t (f vxB-ds. 

Ja Jc Jc 

Now applying Stokes' theorem to the closed contour C: 

[ B-dS = 5t [ V X (v X B) • dS. 

Ja Js 



Hence Ii becomes 



h= [ B-dS- / B-dS = -(5t / Vx (vxB)-dS, 

Js(t+5t) Jsit) Js 



S(t) 

dB 



(2. 



(2.9) 



(2.10) 



IS{t+5t) 

and with ( |2.3|2.6[ ): 

If the induction equation is satisfied, the integrand vanishes and we have shown that 

d^> 



V X (v X B)] • dS 



(2.11) 



dt 



0. 



2.12) 



2.2. Conditions for 'flux freezing' 



The 'freezing' of field lines to the fluid described by (1.27) was a consequence of the 
expression for the electric field, 

E = -vxB/c, (2.13) 
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from which the MHD induction equation followed directly. In deriving this we have made 
the assumption of perfect conductivity. Flux freezing is often described intuitively in a 
rather different way, however, namely in terms of charged particles being bound to the 
field by their orbital motion of around magnetic field lines. Collisions between particles 
would make them jump to orbits on neighboring field lines. Such an 'orbit' picture of 
tight coupling thus applies if the collision rate ojc is much smaller than the cyclotron 
frequency of the charges (also called gyro frequency of Larmor frequency). Such a 
plasma is said to be collisionless or strongly magnetized. 

In this strongly magnetized limit, the electrical conductivity in a direction perpendic- 
ular to the magnetic field would therefore vanish. This is the opposite of the infinite 
conductivity assumed in deriving (2.13). Because of this it is sometimes argued that an 
MHD approximation is not valid in a collisionless plasma. Eq. (2.13) still applies in this 
limit, however, though for a different reason. If an electric field is present and collisions 
negligible, particles drift across the magnetic field with a velocity, independent of their 
charge, of v = — cE x B/i?^ ('E cross B drift'). In a frame comoving with this velocity, 
the electric field vanishes (problem 3.2). If the mean velocities of the different particle 
charged species in the plasma are also representative of the velocity of the fiuid as a 
whole, the observation of an electric field component perpendicular to the magnetic field 
just means we are not observing in the rest frame of the fiuid. Translating to this frame 
the electric field disappears again, just as in the case of perfect conductivity. 

If there is a significant difference in the mean velocity between electrons and ions, 
however, or if the velocity of the charge-neutral component of the plasma relative to 
the charge carriers is important, we are in an intermediate regime between the limits 
of perfect conductivity and perfect magnetization. The MHD induction equation must 



then be extended to include Hall drift and/or ambipolar diffusion terms, see section 2.9 



2.3. Labeling of field lines: Euler potentials 

An explicit labeling of field lines, i.e. numbers that are constant along field lines, is 
sometimes useful. In 3 dimensions such a label requires 2 scalars, counting in two direc- 
tions perpendicular to the field lines. Let xi^it) be a scalar function, taken sufficiently 
smooth, with the property that B • Vx = 0. Isosurfaces x = est. then are called mag- 
netic surfaces: continuous surfaces that are everywhere parallel to the magnetic field. 
The surfaces can be labeled by the value of x- With another scalar ^ define a second, 
intersecting, set of magnetic surfaces, i.e. such that Vx is nowhere parallel to V.^. Each 
field line can then be labeled with its values of x ^-^id ^. 

In general, such a labeling is possible only in a limited region. Field lines on a torus, 
for example, generally wrap around it on 'irrational' paths: without ever returning to 
the same point. To label field lines on such a surface, a cut has to be made across it, 
where the labels are then discontinuous. For any field B there are special scalars Xi i 
such that 

B = Vx X V^; (2.14) 

these are called Euler potentials; they satisfy divB = 0. [Exercise: show this.] They are 
still not quite unique: for a given field B(r) there exists freedom of choice of a single 
(scalar) function along a certain path in space. For more about Euler potentials, see 
Stern (1970). 
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2.4. Reconnection 




Figure 2.2.: Inflow and outflow (black arrows) near a reconnection region (current density in 
green). 

By reconnection we mean the change of topology of a (part of) a field configuration, 
when neighboring field lines running in different directions 'touch' and exchange the 
directions their paths. There are still significant open issues in theories of reconnection, 
hence the subject is not included in chapter 1. 

The simplest and most widely studied case is that of a pressure-supported current 



sheet, as in Fig. 1.16 (upper right panel). It is a model for reconnection in solar flares of 
the classical 'two ribbon' type (cf. Priest 1982 Ch. 10) and the so-called geo-tail of the 
Earth's magnetosphere. The field lines on the two sides of the sheet are antiparallel, and 
the configuration is independent of the y— coordinate (perpendicular to the plane). One 



envisages that finite diffusion allows this sheet to change into a configuration like Fig. 2.2 
There is to be a 'reconnection point' X (actually a line of some finite width, perpendicular 
to the plane) of some extent still to be determined. Because of the symmetry assumed, 
there could in principle be several such points along the sheet, but in less symmetrical 
practice the process will typically start at one point which then preempts reconnection 
elsewhere in its neighborhood. The fluid, with the fleld lines frozen into it, flows towards 
X (top and bottom arrows). Inside X the field lines from above and below are not 
frozen in any more and exchange directions. The field lines newly connected in this 
way are strongly bent. The tension force pulls them away from X (horizontal arrows). 
This causes more fiuid carrying field lines with it to fiow towards X, resulting in a 
quasi-stationary reconnection fiow. 

A logical definition of the rate of reconnection at such a point is the number of field 
lines reconnected, i.e. a magnetic flux per unit time (j) (dimensions G cm^s~^). If the 
extent of the reconnection region in the perpendicular (y-) direction is L, this can be 
written in terms of the velocity Vr at which field lines are advected into the reconnection 



region (vertical arrows in Fig. 2.2): 



= BiLvr, (2.15) 

where is the fleld strength in the inflow region. The task of a reconnection theory is 
then to provide a number for v^- Each of the reconnecting field lines carries an amount 
of fiuid with it that has to change direction by 90° from vertical inflow to horizontal 
outflow. If d is the thickness of the reconnection region, D its extent in the x— direction. 
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and p the density, this amount is D/d ■ p/B, per unit of reconnected magnetic flux. The 
inertia of the mass to be accelerated out of the reconnection region by the tension forces 
thus increases proportionahy to D. The inflow speed and reconnection rate decrease 
accordingly. For more on this subject see Kulsrud Ch. 14.5. 



The sketch in Fig. 2.2 assumes Z) ~ d and a correspondingly high reconnection rate. 
This is in line with indirect evidence such as observations of solar flares and with recon- 
nection in the plasma laboratory. A popular guess adopted in astrophysics is fr ~ 0.1 f a 
where va is the Alfven speed in the inflow region. A problem with this is that such 
a rate is not well supported by the results of numerical MHD simulations that include 
only an Ohmic diffusivity. Instead of forming a small reconnection area, the reconnec- 
tion region tends to spread out into a layer of signiflcant width D S> d (a 'Sweet-Parker' 
configuration) . 

Much better agreement is obtained in MHD simulations of ion-electron plasmas if Hall 



drift (sect. 2.9) is included. If the plasma is sufficiently highly ionized that ion-neutral 
collisions are frequent, ambipolar drift can be ignored. At the high current density in the 
reconnection region the Hall term, however, becomes the dominant term in the induction 



equation (2.26). The induction equation then has the form shown in eq. (2.30), that is, 
the velocity determining the evolution of the magnetic field in the reconnection region is 
that of the electron fiuid. It is effectively decoupled from the ion fiuid in this region. This 
greatly increases the speed of the process, since the field lines can now reconnect without 
having to defiect the heavy ion fiuid through the reconnection region. The coupling of 
the field to the ions becomes important again at some distance from the reconnection 
region. The curvature force of the reconnected field lines drives the ion fiow outward 
there, as in classical reconnection models without a Hall term in the induction equation. 

The reconnection region in an electron-ion plasma becomes effectively 'point-like' 
[d K, D, a 'Petchek' configuration) and fast (infiow at a significant fraction of the Alfven 
speed). For details see Uzdensky &: Kulsrud (2006). This form of reconnection depends 
in an essential way on the Hall effect; in a positron-electron plasma reconnection will be 



different (cf. sect. 2.9). 



2.4.1. Reconnection at low /5 

The situation is intrinsically rather different in a low-/3 environment. A low-/3 field 



is close to a force- free configuration (unless very dynamic, w ~ fA Cg, see 1.4). 



If a current sheet develops in it, it will be closer to the force- free variety than the 



(pressure-supported) current sheet assumed in most reconnection studies (cf. sect. 1.11). 
Since force- free fields are so intimately connected with their boundaries (sect. 1.5.3), 
reconnection in this case is much more strongly determined by the boundaries than in a 
pressure supported sheet. 

A well-studied example of reconnection in force-free fields is magnetic heating of the 
solar corona. The field in the low-density, but highly conducting atmosphere is anchored 
in the high-density plasma of the convective envelope. The fiuid motions there displace 
the footpoints of the magnetic field lines extending into the atmosphere. The displace- 
ments at the two ends of a magnetic loop are uncorrelated, and the field gets wrapped 
into a tangled configuration. Since the displacements are slow compared with the Alfven 
travel time along the field, the wrapping takes place as a series of quasi-static force-free 
configurations. 

Any sufficiently nontrivial displacement process at one or both ends of a bundle of 
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Figure 2.3.: Left (schematic): Field lines of different directions wrapping over (red) and under 
(blue) a force-free current sheet (green). Right: complex current sheet (yellow) 
and field lines wrapping around it in a low-/3 numerical simulation (detail, from 
Galsgaard & Nordlund 1996). 



field lines causes chaotic mixing ( 'http •.//en. wikipedia . org/ wiki/Chaot ic_iiiixing' ) . 
Suppose the imposed displacements at the boundary are smooth in space and in time. 
The direction of the field lines then varies smoothly across the volume. As the wrapping 
process proceeds, however, the length scale on which their direction varies decreases 
exponentially with time. A set of force-free current sheets develops in the volume. When 
the length scale across the sheets reaches the limit where ideal MHD applies, reconnection 
and energy dissipation set in. Because of the exponential nature of chaotic mixing, this 
happens after a finite amount of time which depends only logarithmically on the magnetic 
diffusivity. 

This chaotic mixing process is a kinematic contribution to the rapid development of 
currents sheets (i.e. describable by the induction equation). Concurrent with it there 
is a dynamic contribution, resulting from the Lorentz forces that develop in a tangled 
flux bundle. Locally, the tangling leads to configurations as sketched in the left panel 
of Fig. 2.3 The tension forces in the magnetic field lines wrapping around the two sides 
of a current sheet squeeze it from both sides, much like neighboring links in a chain 
press on each other at their surface of contact. In a force-free sheet, the volume between 
the two directions is filled with field lines of intermediate direction (lower right panel in 
Fig. 1.16). For reconnection to proceed, these have to somehow 'get out of the way' to 
make room for new field lines entering the reconnection region. 

Like all field lines, however, these intermediate field lines are tied to the boundaries. 
This is unlike reconnection in the pressure supported sheet of Fig. 2.2, where the fluid 
leaving the reconnection region does not have such a restriction. Close to a boundary the 
restriction is strong and the formation of small length scales is dominated by the chaotic 
mixing effect of the displacements on the nearby boundary. Further from the bound- 
aries, reconnecting field lines can move out of the way more freely, and current sheets 
become thinner more rapidly (van Ballegooijen 1985). With increasing distance from 
the boundary, this process converges to the 'topological reconnection' picture developed 
by Parker (1979, Ch. 12). 

Numerical simulations of a low-/? field tangled by displacements at a boundary (Gals- 
gaard & Nordlund 1996) show how the magnetic configuration develops a fractal distri- 
bution of current sheets (Fig. 2.3). It takes a finite number of uncorrelated displacements 
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before the sheets are thin enough for dissipation to set in. When these displacements are 
halted in the simulation, dissipation drops quickly by broadening of the current sheets; 
the wrapping pattern of entire field configuration 'freezes'. When the displacements are 
resumed, dissipation immediately springs into action again: reconnection in a low-/3 field 
is driven by the boundaries. 



2.5. Charged clouds 

While ( |1.15[ ) shows that a charge density cannot change in MHD, one might ask what 
happens when a charge density is present in the initial conditions. For example, can one 
still use MHD in a star or a conducting gas cloud that has a net charge? 

In MHD, as in electrostatics, a net charge introduced in a conducting fluid object 
appears only on its surface. If a charge were initially distributed through the star, its 
volume would quickly get polarized by small displacements between the positive and 
negative charge carriers, canceling the charge density in the inside and leaving the net 
charge as a surface layer on its boundary. The time scale for this to happen is governed 
by the plasma frequency, which is usually very fast compared with the MHD time scales 



of interest (see sect. 2.7). Alternatively, this can be seen as a setting a minimum on the 
time scale of the process under consideration, for the MHD approximation to apply. 

Physical charge densities, i.e. net imbalance between plusses and minuses as measured 
in the rest frame of the fluid, thus appear only on boundaries between the conducting 
fluid and non-conducting parts of the system. In frames other than the rest frame of the 
fluid, charge densities appear only as frame transformation quantities, see next section. 



2.6. The charge density in MHD 

In the limit <C c, charge densities in MHD are 'small' in some sense, and are ignored 
because not needed for solving the equations. It is instructive nevertheless to check 
that the procedure followed for deriving the MHD approximation has not introduced an 
inconsistency with Maxwell's equations. In addition, a closer look at charge densities is 
important for interpreting their role in relativistic MHD. 



To test if the third of Maxwell's equations (1.5) is satisfied (which was not used 



in deriving the MHD equations), take the divergence of the electric field under the 



assumption of perfect conductivity, eq. (1.2). Using the formula for the divergence of a 
cross product we have: 

47rcr = V • E = -V • (v X B)/c = B • (V x v)/c- v • (V x B)/c. (2.16) 

(This is valid for arbitrary v/c, since the expression for E is relativistically correct). This 
does not, in general, vanish. Hence we have the disconcerting conclusion that a charge 
density appears in a theory which uses charge neutrality as a starting assumption. 

The key here is the distinction between charge and charge density. To see how charge 
densities appear consider first an irrotational flow (i.e. V x v = 0), so that the first 
term vanishes. The second term describes a charge density that appears in a flow v 
in the direction of the current density j. It can be understood as the consequence of a 
Lorentz transformation between the rest frame of the plasma and the frame in which 



(2.16) is applied. The relativistic current is the four-vector J = (jjCcr). If the primed 



quantities j',<T' are measured in the rest frame of the plasma, and unprimed quantities 
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in the frame in which the plasma velocity is v, a Lorentz transformation yields [with 
/3 = ^/c,7 = (l-/32)-i/2]: 

j =^(^f -i3ca'), (2.17) 
ca =^{ca'-pj'). (2.18) 
If the charge density a' in the plasma frame vanishes, we have j = 'jj' , and 

a = -hj'/c = -I5j/c. (2.19) 



With the current j = cVxB/47r this agrees with the second term in (2.16). 

Thus, while individual charges are Lorentz invariant, the charge density is not because 
it is a quantity per unit of volume, and volumes are subject to Lorentz contraction. In the 
presence of a current, it can be given an arbitrary (positive or negative) value \a\ < j/c 
by a change of reference fram^ 

This formal result can be made more intuitive by taking a microscopic view of the 
plasma, consisting of charge carriers of opposite sign with densities n"*" , n~ , which we 
assume to be equal in the rest frame of the plasma. In the presence of a current, they 
move in opposite directions. A Lorentz transformation with the velocity v then has a 
different effect on the two charge densities, since one is flowing in the same direction as 
V, the other in the opposite. The Lorentz contraction of the two carrier densities differs, 



so they are not the same anymore in the observer's frame (problem 3.28). 



To see what the first term in (2.16) means, consider the plasma in an inertial frame 



locally comoving at some point in space. The electric field, as well as the charge density 



due to the second term in (2.16) then vanish at this point. Unless V x v = 0, however, 
the electric field in the environment of the point does not vanish. It varies approximately 
linearly with distance near this point, implying the presence of a charge density. In 
a frame rotating locally and instantaneously with rotation vector = V x v, this 
charge density vanishes again. By observing in frames rotating at different rates, charge 
densities of either sign appeai[^ 

Charge densities thus disappear in a locally comoving, locally corotating frame of 
reference: the frame in which the plasma is 'most at rest'. An observer in this frame 
sees neither a charge density nor an electric field, and senses no electrical force. Since 
this frame is in general a different one at each point in the fiow, it is not a practically 
useful frame for calculations. It illustrates, however, that charge densities in MHD can 
be regarded as frame transformation quantities. This is the case in general, not only in 
the limit v/c <C 1. Whether relativistic or not, fiows do not create charges. 



What remains to be checked is if the charge density (2.16) needs to be taken into 
account in the non-relativistic equation of motion, in frames other than the fiuid frame, 
since we have ignored this possibility in deriving the MHD equations. The electric field, 
of order v/c, must be kept in order to arrive at the MHD equations. The force density 
~ Ecr is one order in v/c higher, hence can be consistently ignored in all reference frames, 
in the nonrelativistic limit. 

Relativistic MHD is beyond the scope of this text. For a practical formulation and its 
use in numerical applications see Komissarov (1999). 

^ Since — c^a^ is Lorentz invariant, a sufficient condition for this to be possible is ttiat c^a^ < . 
Tills is unlikely to be violated in a medium of any significant conductivity. In the case c^a^ > j^, 
there is a frame in which the current vanishes but a charge density is present. The absence of a 
current in the presence of the electric field implies that the medium is an electric insulator. 

^ Globally rotating reference frames are relativistically problematic. Here we need only a locally coro- 
tating frame in a region of infinitesimal extent, involving arbitrarily small velocities. 
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2.7. Applicability limits of MHD 



A number of different considerations lead to limits of applicability. Start with the condi- 
tion that a charge density is absent (in the frame in which the plasma is 'most at rest', cf. 
sect. 2.6). In the example of a charged cloud (sect. 2.5) this requires that the time scale 
on which a charge density initially present in the volume is canceled by the formation 
of a surface charge, is short compared with the MHD time scales of interest. This time 
scale Te is determined by the most mobile charge carriers (i.e. electrons, usually). In a 
volume of size L, equate the electrostatic energy of the charge density with the kinetic 
energy of the charges, of mass m, when accelerated in the corresponding electric field. 
This yields 

m/(47rne^) = l/w^, (2.20) 

where Wp is called the plasma frequency. This is independent of the size L of the volume. 

In this electrostatic adjustment process the charges that were assumed to be present 
initially do not move to the boundary themselves; this would take longer. Instead, 
their presence is canceled by polarization: small local contractions or expansions of the 
electron fluid relative to the ions. The time of this process is given by 1/wp, or the light 
travel time across the volume, whichever is longer. For the MHD approximation to hold, 
the shortest time scale of interest in the MHD process studied must be longer than this. 



The charge density (2.16) also implies a minimum plasma density that must be present 
for MHD to be applicable: the density of charges implied by the current cannot exceed 
the total density of charged particles n. If the length scale on which B varies in a 



direction perpendicular to itself is L^, the second term in eq. (2.16) thus requires that 
n > vB / [ATvLBec). This is not the most relevant limit involving Lb, however. The 
condition that the current density can be carried by particles moving below the speed 
of light is stricter by a factor c/v. With j = [c/Att) B/L-q = vne : 

n>— V- (2-21) 



If the first term in (2.16) dominates, and is a characteristic value of w/|V x v|, then 
n has to satisfy 

n > vB/{47rL^ec), (2.22) 

or, writing 0, = v/Ly-. 

n>nc, = ^. (2.23) 
This is called the Goldreich- Julian density. If the charged particle density n is pre- 



determined, (2.21) and (2.23) can alternatively be seen as setting minima on the macro- 
scopic length scale Lb and time scale Ly/v for MHD to be applicable. 

The moving charges carrying the electrical current also carry momentum. A change 
of the current in response to a changing field strength requires a force acting on these 
charges. When should one worry about the inertia of the charges responding to this force? 



As a model, take the shearing field configuration of Fig. 1.5 A velocity of amplitude v 
changes direction over a length L perpendicular to it. The current density j, carried 
by particles of mass nic, charge e and number density Uc moving at a velocity Vc is 
j = eUcVc. The momentum pj carried by these charges is thus pj = jnic/e. The current 
in this example increases linearly with time at the shear rate, = v/L. The force 
(per unit volume) needed to achieve this increase is dtPj = {j/Tj){mc/e). To see how 



47 



important this is, compare it with the Lorentz force acting on the fluid as a whole, 
Fl = jB/c. The ratio is: 

^ = ^ = (2.24) 
Fl eBTy ujbTv 

where cob = eB/{mcc) is the cyclotron frequency of the charges, ~ 3 s^^ for 
electrons. The condition that the inertia of a current carried by electrons can be neglected 
is thus 

T„>310"'^-B, (2.25) 
a very mild condition except for time scales deep inside microscopic plasma processes 



(see also sect. 2.8 below). 



Plasma processes of various kind can set in already at densities higher than (2.23), 
leading to an effective increase of the electrical resistivity ('anomalous resistivity'). Even 
when the configuration is well within the limits for ideal MHD to apply on global scales, 
its evolution can lead to the formation of current sheets: regions where the direction of 
the field lines changes so rapidly across magnetic surfaces that field lines are locally not 
'frozen in' anymore. Reconnection of field lines takes place there, with consequences for 
the field configuration as a whole (sects. |2.4[ ). 



Collisionless plasmas 

In many astrophysical fluids the collision frequency oo^ is much lower than the (electron) 
cyclotron frequency ojb'- the plasma is called 'collisionless', or 'strongly magnetized'. 
This differs from the Ohmic conductivity picture that figures somewhat implicitly in 
Chapter 1. This might lead one to conclude that such situations are in the domain 
of plasma physics, not MHD. In the limit lOc ^ however, the charge carriers are 
again tied to the field lines (though for a different reason than under the assumption 



of infinite conductivity, see sect. 2.2). As a result, the MHD equations again apply 
on length scales where local fluid properties like pressure, density and velocity can be 
meaningfully deflned. (Hall drift and ambipolar diffusion terms in the induction equation 



may become relevant, however, see section 2.9) 



For the scales of structures observed in the universe this is almost always the case. 
In the solar wind, for example, the Coulomb interaction length can be of the order of 
an astronomical unit, but the observed structure of its flows and magnetic fields agrees 
closely with expectations from MHD theory. Another example are the structures seen 
in jets and radio lobes of active galactic nuclei, where Coulomb interaction times can 
be of the order of the age of the universe, yet their morphology, with ages of order 10^ 
years, closely resembles structures expected from normal compressible MHD fluids. For 
more on this see Parker (2007, Chs. 1 and 8). 

2.8. The microscopic view of currents in MHD 

With the consistency of MHD established, the subject of currents has become of marginal 
significance ('just take the curl of your field configuration'). This argument does not 
always suffice to dispel concerns how a current can come about without a battery driving 
it. What sets the particles in motion that carry the current? 

The origin of the current of course lies in the fiuid flow. The current in resistive 



MHD results from the flow across field lines (as in Fig. 1.15). In the limit of infinite 
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conductivity, the flow speed associated with this current becomes infinitesimally small, 
but conceptually is still present. Since flow across stationary field lines only occurs at 
finite resistivity, a microscopic view of the source of currents requires a model for the 
resistivity. 

As a simple model which allows for an adjustable resistivity, assume a plasma with a 
small degree of ionization, and consider first the case of high collision frequency (wc ^» 
u)b)- Assume a magnetic field anchored in an object that is stationary in the frame of 
reference. The finite resistivity of the fluid allows it to flow across B with a velocity v 
(cf. sect. 1.10). Between collisions, charges traveling with the flow are displaced from 
the direction of the mean flow by their orbits in the fleld B. Charges of opposite sign 
are displaced in opposite directions, i.e. a current results, perpendicular to v and B. In 
the opposite case of high magnetization, the charges orbit around the field lines with 
infrequent collisions. The collision probability is enhanced during the part of their orbits 
where they are traveling opposite to v. Such collisions reduce the velocity of the orbits 
and at the same time shift their guiding center, again in opposite directions for the +'s 
and — 's. 

For a given velocity, the current produced is proportional to the density of charges 
participating in the collisions, i.e. the degree of ionization in the present example. Con- 
versely, for a given value of the current, the incoming flow speed decreases with increasing 
degree of ionization. The current distorts the magnetic field; the associated magnetic 
force is balanced by the momentum transfer in the collisions. In the limit of vanishing 
resistivity, the fiow speed perpendicular to the field vanishes, in a stationary fleld con- 
flguration. In the general nonstationary case, the magnetic field changes by 'dragging 
with the flow' as described by the MHD induction equation. 

A current needs a flnite time to be set up or change direction. The discussion above 
shows that the relevant time scale is related to the orbital time scale of the charges in 
the magnetic field. See section [2?f| for details. For further discussion of the above in the 



context of a practical application, see section 2.10.3 



2.9. Hall drift and ambipolar difFusion 



The induction equation in the form ( |1.142 ) includes only an 'Ohm's law' resistivity. For it 
to be relevant as the dominant deviation from ideal MHD, differences between the mean 
velocities of neutrals, ions and electrons must be unimportant. When such differences 
do become relevant. Hall drift and/or ambipolar diffusion have to be accounted for. Let 
n, Be, rii be the number densities of neutrals, electrons and ions, and p, p\ the mass 
densities of neutrals and ions, respectively. The MHD induction equation then has the 
form 

5B 

= Vx [(v + VH + Va)xB-r?VxB], (2.26) 

ot 

where v is the velocity of the neutral fiuid component, and 

c V xB j 

VH = Vc-Vi = -- = (2.27) 

47r e ne e 

is the Hall drift: the velocity Ve of the electrons relative to the velocity Vi of the ions. 
The term 

1 (V X B) X B Fl , , 

Va = Vi-V=— ^ ^ = -^, 2.28 

47r 7p/9i 7ppi 
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is the ambipolar drift velocity: the velocity of the ions relative to the neutral component. 
For the derivation of these expressions and for the assumptions involved see Balbus 
(2007). 

The main assumption made is apparent in the form of the expressions on the right 



in (2.27) and (2.28). The second equality in (2.27) shows that the Hall drift is just the 



mean velocity of the electrons associated with the current they carry. This involves the 
approximation that the contribution of the ions to the current can be neglected, which 
requires that the ratio me/m-i of electron to ion mass is neglected. 

The Lorentz force Fl acts directly only on the charged components, since the neutrals 
do not sense the magnetic field. Transmission of the force from the ions to the neutral 
component involves momentum transfer, and associated with it there is a velocity differ- 
ence Va between neutrals and ions. The collision rate is proportional to both the density 
of the neutrals and the density of the ions, with a 'friction coefficient' 7 (with dimensions 
cm^ s~^g~^). This explains the appearance of their densities in the denominator of the 
ambipolar drift velocity (2.28). Both vh and Va vanish at large electron density ng (= ni, 
for singly ionized atoms). 



Eq. (2.26) is the appropriate form of the induction equation when Hall and ambipolar 



drift can be regarded as small flows relative to the neutrals. If the plasma is nearly fully 
ionized, on the other hand, the ions dominate the mass density p. In this case it is more 
useful to write (2.26) in terms of the mean velocity of the ions Vi = v + Va instead of 



the neutrals. With the first equality in (2.28), this yields 



5B 

~dt 



V X [(vi + vh) X B - ??V X B]. 



(2.29) 



The ambipolar process is then just the drift of the inconsequential neutral component 



relative to the ions. With the first equality in (2.27) we can write this also in terms of 
the velocity Vg of the electron fluid alone: 



aB 



V X [vc X B - r^V X B]. 



(2.30) 



This shows that in an ion-electron plasma the magnetic field is strictly speaking not 



'frozen in' in the plasma as a whole, but in its electron component. Eq. (2.30) is not 



a useful form for use in the full MHD problem, however, since the equation of motion 



involves the velocity of the heavy component, not the electrons. Instead, (2.26) or (2.29) 
are needed, depending on the degree of ionization. 

The Hall drift in the above applies to an ion-electron plasma. Its existence depends on 
the mass difference between the positive and negative charge carriers. In a pair plasma^ 
consisting of equal numbers of positrons and electrons, the situation is symmetric, the 
current is carried equally by the positive and negative charges, and the Hall drift vanishes. 
Hall drift probably plays a crucial role in the reconnection of field lines in ion-electron 



plasmas (sect. 2.4); in pure pair plasmas reconnection is likely to function differently. 
This is relevant for the magnetohydrodynamics of pulsar winds as observed in a young 
pulsar like the Crab, which are believed to be consist of nearly pure pair plasmas. 
The ambipolar drift velocity can be decomposed as 



+ Vir, 



(2.311 



where Vgo and Vir are its solenoidal (divvgo = 0) and irrotational (Vxvir = 0) com- 
ponents. The two have different consequences. The irrotational flow component is the 
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'divergent/convergent' one: it causes a gradual pile-up of the charged relative to the 
neutral fluid (as opposed to a rotational motion, which in principle can proceed in a sta- 
tionary state) . Pressure balance of the fluid as a whole implies that a gradient in pressure 
develops between the two, by which Vir eventually stalls, leaving only Vgo (Goldreich and 



Reisenegger 1992). See problem 3.29 



A somewhat implicit assumption in the above is that the density of the charge carriers 
relative to that of the neutrals is determined instantaneously, compared with the MHD 
time scales of interest. This is rarely a concern, since ionization/recombination rates 
are usually very frequent compared with the MHD time scale. It is important, however, 
for the magnetohydrodynamics of a neutron star interior. Relaxation of the charged 
(proton-|-electron) to neutron density ratio by slow weak interaction processes plays a 
role in this case (Goldreich and Reisenegger 1992). 

2.10. Polarization 

2.10.1. Conducting sphere in a vacuum field 



The charge densities discussed in section 2.6 are frame transformation quantities that 



disappear in the locally comoving, corotating frame of the fluid. At the boundary with 
a non-conduction region (a vacuum or an unionized gas) they become real, however, in 
the form of surface charge densities. 

As an example, take a conducting sphere of radius R moving with velocity v into 
an initially uniform vacuum magnetic field B in which the electric field vanishes before 



introduction of the sphere. Assume that the magnetic diffusivity i] (see sect. 1.10) of 
the sphere is small but finite, and observe the sequence of events in a frame comoving 
with the sphere. After its introduction into the magnetic field, the field inside the sphere 
vanishes initially, shielded by currents on its surface. 

The magnetic field diffuses into the sphere on a time scale r = R'^/rj, until for t ^ t 
the field is again uniform everywher^ In the comoving frame the sphere finds itself em- 
bedded in an external electric field E = v x B/c, but inside it the electric field vanishes 
due to the conductivity of the sphere. This is a classical problem in electrostatics. Let 
the magnetic field be in the j;-direction, the velocity v in the y-direction. Polarization 
of the sphere in the electric field displaces the positive and negative charges with respect 
to each other. Inside the sphere the displacement is uniform hence does not cause a 
net charge density, but on the surface of the sphere it results in a charge distribution 
CTg = 3vB cos9/ {4:TtRc) (charge per unit area), where 9 is the polar angle with respect to 
the z-axis. (Exercise: verify this by consulting Jackson E&M Ch. 2.5). 

Outside the sphere, this surface charge creates an electric field equivalent to that 
of a point dipole at the center of the sphere, inside it a uniform electric field which 
just cancels the external field v x B/c. On transforming to the observer's frame, the 
external electric field disappears, and only the field created by the surface charge remains. 
Inside the sphere the electric field then has the expected frame-transformation value 
Ei-nt = —V X B/c. (Exercise: sketch the electric field lines in both frames). 



^ In the laboratory frame we now have a conducting sphere moving across a magnetic field without 
distorting it. This has confused some authors into questioning the whole idea of flux freezing of 
magnetic fields in ideal MHD. Flux freezing and MHD do not apply outside the conductor, however. 
The key is the appearance of a surface charge density. 
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The polarization charge as seen in a comoving frame is the same [up to relativistic 
corrections of order (v/c)'^ due to Lorentz contraction] as in the laboratory frame, but its 
interpretation very different. In the lab frame, the flow v carries two fluids of opposite 
charge across the magnetic field. The Lorentz force is of opposite sign for the two, 
causing them to displace from each other and leaving an excess of +'s on one side and 
— 's on the other. In the comoving frame of the sphere these Lorentz forces are absent. 
In this frame the surface charge is an electrostatic effect, induced on the conducting 
sphere by an external electric field. 



An important application of the above concerns the magnetic field of a rotating neutron 
star, i.e. a pulsar. The surface temperature of a pulsar is so low compared with the virial 
temperature (the temperature needed to cause the atmosphere to inflate to the size of the 
star itself) that near-perfect vacuum holds already some 10 meters above its surface. The 
interior of the star is perfectly conducting in the MHD sense, so that the electric field in a 
frame corotating with the star vanishes in its interior. As in the simpler example above, 
a surface charge appears at the boundary with the vacuum. In the corotating frame, 
it can be seen as polarization resulting from the external electric field (fi x r) x B/c in 
which the conducting star finds itself embedded. In the inertial frame, where this field 
does not exist, the polarization can be seen as due to the differential Lorentz force acting 
on the two charged species rotating in the magnetic field of the star. 

The surface charge develops at the level in the atmosphere where vacuum conditions 
take over from the high-conductivity regime in deeper layers. This happens at such a 
low density that the magnetic field completely dominates over the gas pressure there, 
f3 <^ 1. Unaffected by the fluid pressure, the magnetic held in the surface charge region 
is thus the continuation of the vacuum fleld outside. 

Suppose the magnetic field of the star is stationary in an inertial frame, 9B/(9t = 0. 
It follows that it must be axisymmetric with respect to the axis of rotation, and by 
the induction equation the electric field has a potential, E = —'Vcj). It therefore has a 
potential everywhere, in the conducting region below the surface charge as well as in 
the surface charge region itself and in the vacuum outside. The radial component of 
E makes a jump across the surface charge, but its potential is continuous. To find the 
electric field in the vacuum, solve V^(/> = there, with the boundary condition that (p is 
continuous across the surface charge. 

As an example take for the magnetic field configuration a dipole aligned with the 
rotation axis and centered on the star. If Bq is its field strength on the rotational north 
pole of the star, the field is 



where r, 6 are the spherical radius and polar angle, z and f unit vectors along the rotation 
axis and in the (spherical) radial direction respectively, and R the radius of the star. 
The electric field Ei in the star, as seen in an inertial frame: 



2.10.2. Pulsars 




(2.32) 




(2.33) 



becomes, with (2.32) 



Ei = -Bo 




i2 6'-l)f-2cos6'z]. 



(2.34) 
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Its potential is (exercise: show this) 

Q O E> 

0i = 5oi?— -(cos2 0-l) + A;, (2.35) 
2c r 

where k is an arbitrary constant. The angular dependence of the field is that of a 
quadrupole. Since the potential is continuous across the surface charge at r = i?, the ex- 
ternal field has the same angular dependence (but a different dependence on r). Matching 
sX r = R yields, (taking the potential to vanish at infinity): 

(^e = i?oi2^(-)'(cos2 0-1/3). (2.36) 
2c r 

The surface charge density cTs is then 

as = (Ee-Ei)-f/47r = -9^(0e-'/'i)/4^ (2.37) 
= --^ cos^e. (2.38) 

(Goldreich and Julian 1969). For the typical magnetic field strengths and rotation rates 
of pulsars the electric field is so strong that the external volume gets filled with charges 
pulled from the surface of the neutron star. This leads to much more difficult problems, 
the subject of pulsar magneto spheres. 

2.10.3. Electricity from MHD 

Experimental devices for the direct extraction of electric power from a gas flow (MHD 



power generation) provide a nice illustration of the discussions in sect. 2.8 and the 
sections above. A hot gas, made partially conducting by seeding it with atoms of low 
ionization potential like Potassium, is made to flow across an externally applied magnetic 
field. On the sides of the flow channel there are electrodes, parallel to both B and v. 
With the lab environment being non-conducting, the setup is similar to the polarization 



experiment in 2.10.1 If the electrodes are not connected by an external circuit, an 
electric field — v x B/c is present in the lab frame, and a corresponding polarization 
charge appears on the electrodes. The current discussed in |2.8| is absent. The defiection 
of the charges in their orbits around the field that would drive this current is canceled 
exactly by the E x B drift due to the electric field in the lab frame. In the fiuid frame, 
the electric field vanishes. The fiuid is not aware that in the lab frame it is seen as 
flowing across a i?-field. 

We now change conditions by closing the external circuit with a negligible resistance. 
The polarization charge is shorted out, and with it the electric fleld in the lab frame 
disappears. In the fluid frame there is now an electric field -|-v x B/c, and associated 
with it is a current. In the lab frame, the source of this current is as described in 2.8 in 
the fiuid frame on the other hand it is just a current driven by an electric field, across the 
internal resistance of the plasma. The flow pattern and the distortion of the fleld lines 



follow from the MHD equations of motion (1.18) and induction (1.142). The flow does 



work against the curvature force of the distorted field lines. In a short-circuited state 
like this, the associated power is all dissipated in the flow itself (like in a shortcircuited 
battery) . 

As this discussion illustrates, it is important to be aware how much physical descrip- 
tions in terms of microscopic processes can differ in the lab and in the fluid frame. One 
has to keep them clearly separate to avoid confusion. 
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2.10.4. Critical ionization velocity 



An interesting situation arises if the moving conductor of section 2.10.1 is replaced by an 
insulator in the form of an unionized (atomic or molecular) gas. Assume again that no 
electric field is present in the laboratory frame. Let the gas consist of atoms A of mass 
m flowing with speed v across the magnetic field. Now insert an ion A+ into the gas 
by hand. Being charged, it is tied to the magnetic field, and its guiding center initially 
at rest in the lab frame. The neutral gas flowing past it collides with the ion, with a 
center-of-mass energy of mv'^/2. The collision can lead to the further ionization of an 
atom of the incoming gas if the collision energy exceeds the ionization energy e of A. 
The new ion will also be tied to the magnetic field, so there are then additional ions 
colliding with the neutral gas, each producing more ions. Runaway ionization of the gas 
will therefore result when v sufficiently exceeds the critical ionization velocity Vc, 

Vc = (2e/m)^/2_ (2.39) 

This yields Vc ~ 40 km/s for A=Hydrogen and ~ 5 km/s for Potassium (one of the 
first elements to contribute a significant degree of ionization in a gas of standard cosmic 
composition). Since the presence of a stray ion is usually plausible, the conclusion is 
that a neutral gas moving across a magnetic field will ionize when its velocity (in the 
frame where E = 0) exceeds Vc (a suggestion due to Alfven). 

An assumption in the simple picture above is that the ions can be considered tied to 
the field, in the laboratory frame, while the neutrals flow past it. This means that the 
density must be low enough for ambipolar drift to be significant: the ambipolar drift 



velocity fa must exceed (2.39). Since the value of v^, decreases as the density of ions 
increases, the runaway is limited: the degree of ionization will eventually saturate at a 
level where Va, is of the order Vc- 



2.11. Surface stress: example 

As an example for evaluating forces by using the stress tensor, consider support against 
gravity by a magnetic field [Fig. |2.4[ a model for the (partial) support of an accretion 
disk by a magnetic field embedded in it]. The acceleration of gravity acts towards the 
left, g = —gx. Mass is present in a layer of thickness 2h around the midplane z = 0. 
Let L = |B|/f?x.|B| be a characteristic length on which B varies along the plane, and use 
a unit of length such that L ~ 0{1). We take the layer to be thin, /i <C 1 in this unit. 
Outside the layer {\z\ > h) the field is assumed independent of the y-coordinate, and 
symmetric about the midplane, Bz{—h) = Bz{h), Bx{—h) = —Bx{h). Inside the layer 
the field can be of arbitrary shape. 

The question to answer is how large the field strength must be to support the disk 
significantly against gravity. We do this by calculating the magnetic stresses over a fiat 
box of unit surface area, — /i<2;</i, 0<x<l,0<y<l. The magnetic force on 



the mass in this box is then given by the right hand side of (1.47). It has contributions 
from the top, bottom, and side surfaces, and the outward normal is to be taken on these 
surfaces. The top and bottom surfaces have normals n = (0, 0, 1), (0, 0, —1), respectively 
and yield contributions 

Fx,top = Mx,{h) = Bx{h)B,{h)/4Tr, (2.40) 
Fx,hot = -MxA-h) = -Bx{-h)B,{-h)/47r = Bx{h)B,{h)/4TT, (2.41) 
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Figure 2.4.: A thin accretion disk with magnetic support against the gravity of a central mass 
(located to the left). 



The contributions from the sides are smaller by a factor h on account of the assumed 
thinness of the layer, hence will be ignored. Summing up, the magnetic force on the 
disk, per unit surface area in the x— direction is 

F^ = iB,Bj2Tr (2.42) 

(where is evaluated at the top surface of the layer). The force of gravity Fg is found 
from the mass in the box. Let 

rh 

p dz (2.43) 



i-h 

be the surface mass density of the layer. The force per unit surface area of the disk is 
then 

Fg = -x^S. (2.44) 

For the magnetic field to contribute significantly to support (compared with the rota- 
tional support due orbital motion), we must have 

Ba:B^ ~ (2.45) 

implying also that B^B^ > {z = h), i.e. the field is bent outward, as sketched in 
Fig.[23 

The example shows how by evaluating the stress on its surface, we have been able 
to circumvent the details of the magnetic forces inside the box. Of course this gives a 
useful result only if the magnetic fields on the top and bottom surfaces are known. In 
the applications we have in mind, vacuum conditions hold outside the layer, and B is a 



potential field there, V x B = (sect. 1.5.2). It is determined only by the distribution of 



field lines crossing the disk, i.e. the vertical component Bz{x,z = h). With B^ given as 
a function of x, the magnetic potential 4){x, z) above and below the disk can be found, 



and from this, Bx = —dx<p and the force (2.42) 



Rotating the figure by —90°, it can also be seen as the sketch for a quiescent promi- 
nence: a sheet of cool gas floating in the solar atmosphere, supported against gravity by 
a magnetic field (see Ch. 11 in Priest 1982). 

2.12. 'Compressibility' of a magnetic field 



Suppose the volume of the box of uniform field in Fig. 1.6 is changed by expansions in 



X, y, and z, and we ask the question how this changes the magnetic pressure in it. Let 
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the expansion take place uniformly, i.e. by displacements Ax ~ x, Ay ~ y, Az ~ z. 



The magnetic field then changes in magnitude, but remains uniform (problem 3.18). 
Consider next an expansion in x and y only, i.e. Az = 0. Under such a change the 
volume V of the box is proportional to its cross sectional area A in the (x, y)-plane. By 
conservation of magnetic flux, the number of field lines crossing A remains unchanged, 
hence i? ~ 1/A ^ 1/^- The magnetic pressure thus varies as B'^/Sir ~ V'"^. 

On adiabatic expansion, the pressure of a gas with a ratio of specific heats Cp/cv = 7 
varies as p ~ V~'^. We can therefore say that, for expansions perpendicular to the field, 
the magnetic pressure varies like that of a gas with 7 = 2; somewhat stiffer than a fully 
ionized ideal gas (7 = 5/3). The tension in the field, however, makes the analogy with 
a compressible gas inapplicable for expansions along B. 

Consider next the case when the box contains a magnetic field that is of unspecified 
configuration, except it is known or arbitrarily assumed that the 3 components are of 
the same average mean square strength. That is, if <> indicates an average over the 
volume, 

< Bl >=< Bl >=< Bl > . (2.46) 

Expand the volume uniformly in the x-direction by a small amount, such that the new 
position xi of everything in the box is related to its old position xq by xi = (1 + e)x. 
Since the induction equation is linear in B, the effect of the expansion can be done 
separately on the 3 components of B. Bx remains unchanged, Bix{xi) = Box{xo), but 
the y- and z-components are expanded in the perpendicular direction, and are reduced 
by the factor 1 + e. To order e the magnetic pressure thus changes by the factor 

Bf{xi,y, z)/Bi{xo, y, z) = [1 + 2/(1 + ef]/3 = 1 - ^e- (2.47) 

The same factor applies for small expansions simultaneously in several directions (prob- 
lem 3.18 i). When the expansion is isotropic and uniform in all directions, ( 2.47| ) even 



holds for fields of different strength in the 3 directions (problem 3.18| 3). 



This is the basis for the statement sometimes encountered that a 'tangled' magnetic 
field behaves like fluid with 7 = 4/3, analogous to the pressure of a photon field. The 
statement is nevertheless of very limited use. Volume changes in a fluid are rarely 
even remotely isotropic (in the incompressible limit: never). Where they are not, the 
strengths of the fleld components rapidly become unequal on distortions of the volume, 



after which the factor (2.47) does not apply anymore. Examples were given in section 
03 
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Figure 2.5.: Jet from a 3D MHD simulation by Moll (2009). Rotation applied near the origin 
(inset) of an initially radial magnetic field (thin lines) produces an outflow with an 
increasingly twisted field that starts buckling under the pressure of the azimuthal 
field component. Ahead and surrounding the jet a fast mode shock (dashed) prop- 
agates into the surrounding magnetic field. The volume current in the interior of 
the jet is compensated by the surface current along the interface with the surround- 
ing magnetic field (green arrows). The pressure of the azimuthal field component 
causes the jet to expand into the surroundings. 
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2.13. Twisted magnetic fields: jets 



An application where twisted magnetic fields play an essential role is the production of 
fast outflows from magnetic fields anchored in rotating objects (stars, accretion disks, 
black holes). 

A rotating disk which contains an ordered magnetic field accelerates an outflow consist- 
ing of a highly twisted field (Fig. 2.5 ). The pressure B'^/Stt contributed by the azimuthal 
(twist) component of this field causes it to expand into the pre-existing field configuration 
(cf. sect. 1.6.1), and produces a shock wave around and ahead of the jet. Near the base 
of the jet the twisted configuration is stabilized by the surrounding field. As the twist 
angle increases with distance from the source, and the pressure of the external medium 
declines, the field configuration of the jet starts buckling under the pressure of the az- 
imuthal component. Since the poloidal field strength (~ 1/t^) declines more rapidly 
with distance r than the azimuthal component (~ 1/r), this happens inevitably at 
some finite distance from the source. 

The inset in Fig. 2.5 shows how the 'closing of currents' in a jet is not an issue. The 
current is not needed for understanding a magnetorotationally powered jet in the first 
place, but if one insists it can be computed from the field configuration. It closes through 
the head of the jet and along its boundary with the environment (compare with section 



1.6.1 and problems 3.15 3.16). 



The interest in magnetically powered jets has given rise to a thriving subculture in 
the literature fueled by misunderstanding of the role of currents. Currents are confused 
with the jet itself, the direction of the currents confused with the direction of the twist 
in the magnetic field, and unnecessary elaborations made on the closing of currents. In 
this tradition currents are correctly understood as due to the rotation at the base of a 
jet, but are then incorrectly dissociated from this process and reinterpreted as the source 
of the toroidal field, thereby mixing up cause and effect. These misunderstandings can 
be avoided with a grasp of the basics of MHD. 



2.14. Magnetic helicity and reconnection 

Since a field of zero energy also has zero magnetic helicity, a field configuration with a 
finite helicity cannot decay completely, however far out of equilibrium or unstable. At 
least not as long as flux freezing holds, since helicity is conserved only in ideal MHD. 
In practice, relaxation of the conflguration can quickly lead to the formation of length 
scales small enough for reconnection to become effective (current sheets), after which 
energy and helicity can both decline (at a slower rate). 

Reconnection, however, does not necessarily cause the helicity of a configuration to 
decrease. Helicity can also be created (quasi out of nothing) by a change of topology. 



Fig. 2.6 shows an example of two disconnected planar magnetic loops, embedded in a 
field- free environment. They have the same amount magnetic flux, both are untwisted 
hence have zero helicity. Assume they are brought together by a fluid flow until they 
touch. (This can be achieved with negligible input of energy). Reconnection between 
the loops releases magnetic energy, but the resulting single loop now has a finite helicity. 
It is contained in the shape of its path (writhe). Folding this loop out onto a plane, it 
is seen that this amount of helicity corresponds to a twist of one turn around the plane 
loop. [The unfolding step requires input of energy, but does not change helicity]. 
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Figure 2.6.: Helicity generation by reconnection. Two untwisted loops (a) with zero helicity 
reconnect at the point of closest approach (b). The resulting single loop is helical; 
on folding it out into a plane the helicity appears as a twist of one turn (c). 



The inverse sequence c-b-a in Fig. 2.6 is of course also possible. One can imagine a 
twisted loop as in c to become unstable, untwisting itself by folding into a figure-eiglit, 
which subsequently reconnects into two untwisted 0-rings. The sequence in Fig. |2.6| is 
more likely in an environment where fluid flows are sufficiently dominant {v/va ^ 1) 
to bring initially disconnected loops into contact. The inverse could take place in a 
hydrodynamically more quiescent environment {v/vx ^ 1), where the sequence of events 
is dictated instead by the magnetic instability of the twisted loop itself. 

The field configuration of Fig. 2.6 can be enclosed entirely in a simply connected 
volume V , so its helicity can be defined uniquely (and is conserved until reconnection 
takes place). Intuitively one feels the need for something that can be defined also for 
cases where B • n does not necessarily vanish everywhere on the surface of the volume we 
are interested in. With some restriction, this is possible in terms of a relative helicity. 

If the normal component B • n does not vanish, but can be kept fixed on the surface 
of y, we can complement the region outside V by attaching to it an external volume 
Vq. In this volume we imagine a fictitious magnetic field Bg that 'connects the ends' of 
the field lines sticking out out V . That is: a) on the part of the surface shared by V 
and Ve the normal component Be • n matches that of B on the surface of V (so that 
divB = is satisfied everywhere, b) Be • n = on the remaining surface of 14, and c) 
ideal MHD holds inside Vf,. With this artifice, the helicity is still not defined uniquely 
since it depends on the precise choice of Bg, but changes in helicity are now uniquely 
defined irrespective of the shape of Be (Berger &: Field 1986). That is, the helicity of 
the entire volume can now be meaningfully attributed, apart from a fixed constant, to 
the helicity of B in V only. Essential for this to work is that B • n does not change on 
the interface between V and Vf, and that this interface stays fixed, whatever happens 
inside V or is imagined to happen in Ve- 

This is illustrated with the example sketched in Fig. 2.7, a variation on the inverse 
sequence in Fig. 2.6 It shows a simplified view of events that can happen on the surfaces 
of magnetically active stars like the Sun, or strongly magnetic neutron stars (magnetars). 
A section of a twisted magnetic flux bundle initially buried below the surface rises 
through the surface (by magnetic buoyancy for example, cf. section 1.3.5). Only the 
erupting segment interests us, but it has open ends. In order to measure the changes in 
helicity in this sequence, we attach a fictitious bundle of field lines to it (a 'spectator field', 
dashed lines in the figure) so as to complete the segment into a loop. This configuration 
can then be embedded in a volume where helicity is defined uniquely. 

The low density environment of the atmosphere of the star does not support twist to 
the same degree as the high-density interior. The field reduces twist by exchanging it 
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Figure 2.7.: Emergence of a twisted field through the surface of a star. Left: a section of a 
strand below the surface (black) contains one turn of twist. Middle: it erupts 
through the surface, while its path makes a half-turn rotation. Right: an untwisted 
loop reconnects from it into the atmosphere while leaving a configuration of zero 
helicity. Dashed: loop completion used for the definition of helicity. 



for a loop in its path. Helicity is conserved in this step. By reconnection at the point 
where the path touches itself (middle panel) a loop of field disconnects, and is ejected 
from the atmosphere by the pressure gradient in the surrounding magnetic field (see 
problem 3.20). The magnetic bundle is now contained below the surface as before, but 
is untwisted, and the helicity of the configuration has decreased by the equivalent of one 
turn of twist. 

The helicity has disappeared already at the moment of reconnection, while nothing 
helical has left the star. This is an idealized example, but it serves to illustrate that the 
change of helicity in a magnetic ejection process like Fig. 2.6 is not properly described 
by a term like 'helicity ejection'. Current helicity (sect. 1.6.2) does not fare much better. 
Unlike magnetic helicity, this is a quantity that does have a local definition, but since 
there is no conservation equation for it like there is for mass or energy, it does not make 
sense to describe the ejection process in terms of a 'flux' of current helicity either. All 
one can say is that unwinding of twist that was present in an emerging loop of field 
lines leads to the formation of a plasmoid that is ejected from the star, while helicity 
decreases in the process. In the Sun, such processes are responsible for the coronal mass 
ejections that cause aurorae and disturbances in the Earth's magnetic field. 

For a more on the properties of magnetic helicity, and its relation to vorticity in 
ordinary fluid mechanics see Moffatt (1985). 
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3. Exercises and problems 



3.1. Currents from flows 

This exercise gives an example how currents appear and disappear when a flow acts 
on a magnetic field. An initially uniform magnetic field B = i?z in the z-direction is 
embedded in a flow in the x-direction, varying with z as v = tanh(2;) x. From the MHD 
induction equation calculate how the magnetic field has changed after a time t. Calculate 
the current distribution. Now reverse the sign of v and notice how the currents have 
vanished again after another time interval t. 



3.2. Particle orbits 

Calculate the orbit of a charged particle, initially at rest, in a magnetic field with B in 
the y-direction and E in the z-direction. Show that the motion consists of a circular 
motion around the y-axis superposed on a uniform velocity Vd = cE/B in the x-direction. 
Sketch the resulting orbit. The velocity is called the E x B drift velocity. Note that 
it is identical to the fiuid velocity v in (1.2), hence the same for particles of all mass and 
charge. 



3.3. Displacement current at finite conductivity 



In section 1.1.5 perfect conduction was assumed. If an Ohm's law conductivity (1.140) 



is assumed instead, verify with (1.141) that Maxwell's equation (1.3) can be written as 



4^j + |^(-) 

Ot (Jc 



d_ 
di 



V X B/c) = cV X B, 



(3.1) 



in the nonrelativistic limit. As in sect. 1.1.5 the third term on the left is negligible 



compared with the right hand side. With characteristic values for velocity V, length L 



and time scale L/V as in section 1.10, show that the LHS in (3.1) is of the order 



4^j(l 



c Lc' 



(3.2) 



where rj is magnetic diffusivity. If I and v are the mean free path and velocity of the 
current carriers, 7] ~ Iv. For the MHD approximation to hold, the macroscopic length L 
must be ^ I. The Ohmic term in the displacement current is thus of order {Vv/c^){l/ L), 
and can also be neglected. 

3.4. Alternative form of the induction equation 



a. Derive eq. (1.25) by using the induction equation in the form (1.24) and the continuity 
equation. 
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b. A horizontal magnetic field B = iJo x of 10^ G lies at the base of the solar convection 
zone, where the density is p = 0.2 g/cm^. If this field were to rise to the surface of the 
Sun {p = 210^^), and assuming it remains horizontal in the process, what would its 
field strength be? 



3.5. Integrated induction equation 

Derive eq. ( |1.75[ ) by integrating the continuity equation from to 5t, keeping only first 
order in small quantities. 



3.6. Stretching of a thin flux tube 

A straight bundle of field lines with strength B{t) oriented along the x-axis is embedded 
in a field- free medium with constant pressure pe- A section with initial length Lq = 1 
of the bundle is stretched uniformly in the x— direction to a new length L{t) = at, 
with a = est. During the stretching, the bundle is in pressure equilibrium with its 
surroundings, pi + B'^/Stt = p^, where pi is the internal gas pressure, p\ = TZpT. The 
temperature T is kept constant in time. Using conservation of mass in the lengthening 
segment, calculate B{t), p{t) and the cross section of the bundle, starting from initial 
conditions po, Bq such that Bq/Stt <^ P^. Describe the initial {t — >• 0) and asymptotic 
{t — )• oo) dependences of B, p. What is the value of the plasma-^ at the transition 
between these regimes? 



3.7. Magnetic flux 

Using the divergence theorem, show that the magnetic flux passing through a surface S 



bounded by a given fluid loop, as in Fig. 1.2 is independent of the choice of S. 



3.8. Magnetic forces in a monopole field 

This exercise illustrates the limitations of identifying the two terms in eq. ( |1.37| ) with 
respectively a pressure gradient and a curvature force. In spherical coordinates (r, 9, cp), 
the held of a monopole is (l/r^,0,0). 

a. Show that the Lorentz force in this held vanishes except at the origin. 



b. Calculate the 'pressure gradient' and 'curvature' terms of eq. (1.37). 

[A magnetic monopole does not exist, of course, but a physically realizable field can be 
made by reversing the sign of the field in the southern hemisphere {6 > tt/2): the 'split 
monopole' configuration. Since the Lorentz force does not depend on the sign of B, this 
makes no difference for the calculation.! 



3.9. Magnetic forces in an azimuthal magnetic field 

In cylindrical coordinates {w, (j), z), consider a purely azimuthal magnetic field B = B (f), 
such that B = Bq/uj, where Bq is a constant. Calculate the magnetic pressure gradient, 
the curvature force, the Lorentz force, and the current along the axis. 
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3.10. The surface force at a change in direction of B 



A field B = i?o z changes direction by an angle a at a plane surface z = zq. Calculate 
the force vector exerted on the fluid on this surface. (Divide the surface into two parallel 
ones separated by an infinitesimal distance over which the change in direction occurs. 
Such surface forces occur in MHD shocks). 



3.11. Magnetic energy and stress 

A spherical star has a uniform magnetic field B in it. Outside the star, the field continues 
as a potential field. 

a. Remind yourself that the outside field is that of a point-dipole centered on the star (cf. 
Jackson E&M, Ch. 5.10). 

b. Calculate the magnetic stress acting at the surface of the star. 

c. Calculate the total magnetic energy of i) the field inside the star, ii) the external field. 

d. A field configuration like this in a star is unstable. Estimate the growth time scale for 
this instability in a star of mean density p (Flowers & Ruderman 1977). 



3.12. Expanding field loop in constant density fluid 



[Exercise in vector calculus, see ( |4.15[ )-(4.25 )]. A circular loop of field lines of radius R 
lies in a perfectly conducting fluid of constant density. It is centered on the z— axis in 
the z = plane of cylindrical coordinates {w,(p,z). The radius expands by an axially 



symmetric flow in the -oj-direction carrying the loop with it. Using (1.25), show that the 
field strength of the loop increases as B ^ R. 



3.13. Magnetic buoyancy 

A horizontal bundle of magnetic field of diameter d and field strength B, in temperature 
equilibrium with its surroundings, rises by magnetic buoyancy in a plasma with density 
p and pressure p. Calculate the speed of rise of the tube in the presence of hydrodynamic 
drag with drag coefficient Cd (look this up in wikipedia). Express the speed in terms of 
an Alfven speed and the pressure scale height H = p/{pg)- (In this exercise, ignore the 
effect of a stabilizing density gradient. See the next problem for the effect of a stable 
stratification). 



3.14. Speed of buoyant rise 

(Exercise in astrophysical order-of- magnitude estimates). 

Consider a magnetic field in a stably stratified star (a main sequence A star, say). As in 
the physics of the Earth's atmosphere, the stability of the stratification is characterized 
by the difference between the temperature gradient and an adiabatic gradient. If p{z) is 
the (hydrostatic) distribution of pressure with depth z (counted positive in the direction 
of gravity) these gradients are, in astrophysical notation (see Kippenhahn et al. 2012 for 
details) : 

^ dlnT ^ ,dlnr (3 3) 

dlnp ' dlnp ^ ' 
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where Va is the dependence on pressure of the potential temperature: the temperature 
after adiabatic compression/expansion from a fixed reference pressure to the pressure p 
of the stratification. It is a function of the thermodynamic state of the gas only. In a 
stable stratification, Va > V. A fluid element displaced from its equilibrium position in 
the stratification oscillates about it with the buoyancy frequency N, given by 

A^' = |(Va-V), (3.4) 

where g the acceleration of gravity, and H = (dlnp/dz)~^ the pressure scale height. If 
it is displaced (in pressure equilibrium with its surroundings) by a small amount 6z, it 
develops a density difference 

^ = ->.-V). (3.5) 

In this stratification imagine a horizontal flux strand of strength B and radius r, initially 
in temperature equilibrium with its surroundings, so it is buoyant (section [1.3.5 ). 



a. From (1.56) calculate the displacement Sz of the strand after it has settled to a new 
equilibrium position, and calculate the temperature difference it has developed with its 
surroundings, assuming the displacement took place adiabatically. Now let the temper- 
ature difference equalize again. The time scale for this, the thermal diffusion time, is 
r = r^/^T, where kt is the thermal diffusivity of the gas. 

b. Estimate the speed of rise of the tube when thermal diffusion is taken into account 
(explain why this is a good estimate). 

c. Express the rise time (tr) over a distance of order of the radius R of the star in terms 
of its thermal relaxation time, the Kelvin-Helmholtz time txH = R'^/kt- 

d. A typical value of txH for this A star is 10^ yr, its central pressure is of the order 10^^ 
erg/cm^. Compare tr with the main sequence life time of the star (~ 10^ yr), assuming 
the largest field strength observed in an A star, ~ 10^ G, and a length scale r/H = 0.1 
for the field. Conclusion? 



3.15. Pressure in a twisted flux tube 



Take a fiux tube like in Fig. 1.8 with an initially uniform vertical field Bq: 

Bo = Soz, (^n<l), Bo = 0, (ro>l), (3.6) 

where w is the cylindrical radial coordinate. We twist it by applying uniform rotation 
in opposite directions at top and bottom, so the velocity applied is 

v = $7zx-c^ (z = l), v = — JIzx-zat {z = —1). (3.7) 

Assume = outside the cylinder and constant inside. To prevent the tube from 
expanding due to the azimuthal field component that develops, we add a compensating 
gas pressure Ap to keep the configuration in pressure equilibrium with Bz unchanged. 
Calculate Ap{w,t). 
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3.16. Currents in a twisted flux tube 



Starting with a flux tube as in 3.15, apply a twist over an angle 6, which depends on the 
distance from the axis: 

9 = 00 (ro<0.5), e = eo{l-uj) (0.5<ro<l), 6^ = (ro > 1). (3.8) 

Assume 1 so radial expansion can be neglected. Sketch, as functions of zu: the 

azimuthal field B^, the vertical current jz and the total current contained within a circle 
of radius zu (compute these quantities if you are not sure about the sketch). Note that 
the direction of the current changes across the tube while the direction of twist does not. 



3.17. Magnetic stars 

[From the literature on Ap stars] Consider a star with a static magnetic field (a magnetic 
Ap star, for example). Outside the star there is vacuum, the field there a potential field 
V X B = 0. Some force-free field is constructed in the interior of the star. The normal 
component of this force-free field is evaluated on the surface of the star. This used as 
boundary condition to construct the field outside the star, using standard potential field 
theory. What's wrong with such a model for Ap stars; why is it not a counterexample 
to the vanishing force-free field theorem? 



3.18. Magnetic compressibility 

(From section [2. 12[ ). 

a. Show that under a uniform expansion, a uniform magnetic field stays uniform. 

b. Next consider a uniform, isotropic expansion, such that the positions of a fluid particle 
initially at ro becomes ri = Ktq. Show that the magnetic pressure varies like K~'^/^, for 
arbitrary initial field Bo(ro). [Be warned, however, of the limited practical significance 
of this fact.] 



3.19. Winding-up of field lines in a differentially rotating star 

A spherical fiuid star of radius R rotates with an angular velocity r2(r) which is constant 
on spherical surfaces but varies with distance r from the center, such that 

Q{r) = no{2-r/R), (3.9) 

At time t = the star has in its interior r < R a uniform magnetic field Bq = Bqz along 
the rotation axis. Outside the star is a vacuum magnetic field. 

a. Remind yourself (or verify with Jackson Ch. 5) that the field outside the star that 
matches to the uniform interior field is a point-dipole of strength /i = BqR^/2. 

b. Using cylindrical coordinates {w, (p, z) with axis z taken along the rotation axis, calculate 
how the magnetic field changes with time. [If you feel like it, also calculate the current 
in this field.] 

c. Note that the external vacuum field does not change since the normal component of B 
does not change on the surface of the star. [If in b you did, why does the current not 
change the field outside the star?] 



66 



3.20. Diamagnetic forces 



An unmagnetized, perfectly conducting sphere is placed in an initially uniform external 
field Bo in the z— direction. 

a. Calculate the change in the external field after introduction of the sphere, assuming the 



field remains current-free (hint: compare with problem 3.11). Sketch the field lines. 

b. ('squeezing' by magnetic forces). Calculate the forces acting on the surface of the sphere. 
If the sphere consists of a fluid of fixed density p, estimate (by order of magntitude) the 
time scale for the sphere to change shape under the influence of these forces. 

c. ('melon seed effect') Replace the initially uniform external field of a) by a potential 
field that decreases in strength along the axis z on a characteristic length scale L. Show 
(qualitatively) that the sphere now experiences a net force acting on it (in which direc- 
tion?) and estimate the acceleration (of the center of mass) under this force. Inclusions 
of low field strength in the fluid thus behave like diamagnetic material. 



3.21. Helicity of linked loops 

Using the conservation of magnetic helicity, show that the helicity of the linked loop 
configuration in Fig. 1.9 is 2<I>i<I>2. (hint: deform the loops, without cutting through 
field lines, to give them shapes for which the calculation becomes easy). 



3.22. Stream function in a plane 

Analogous to the axially symmetric case, define a stream function for a magnetic field 
with planar symmetry (e.g. in Cartesian coordinates x, y, z with dy = 0), and show that 
in this case it is equivalent to a vector potential of the field. 



3.23. General form of a plane Alfven wave 

Show that 



Vy/VA 
dBy/B 



fi{x)gi{z + VAt) + h{x)g2{z 
fi{x)gi{z + VAt) - f2{x)g2{z 



VAt), 
VAt), 



(3.10) 



where /i, /2, 51,52 are arbitrary functions, is a solution of eqs. ( |1.107 ) and (1.109). Note 
that /, g must be continuous, but their first derivative need not be (the field lines can 



have a sharp bend, as in Fig. 1.10) 



3.24. Torsional Alfven waves 

a. Write the linearized equations of motion and induction in cylindrical coordinates {uj, (p, z), 
with B along the z-axis as before. Show that there exist 'cylindrical' Alfven waves, with 
all perturbations vanishing except 5B^ and Vfj,. These are called torsional Alfven waves. 

b. In a perfectly conducting fluid of constant density p there is at time t = a uniform 
magnetic field parallel to the axis of a cylindrical coordinate frame, B(ti7, cj), z,t = 0) = 
Bq z. The field is 'anchored' at z = in a plate of radius R in the (ro, (j)) plane. At 
t = it starts to rotate uniformly around the z-axis with rate VL. Assuming this can 
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be treated as a linear perturbation, what is the field line angle in the torsional 

Alfven wave launched by this setup? Calculate the torque T{w) acting on the plate, 
and the work done by the plate against this torque. 

3.25. Currents in an Alfven wave 

Sketch the electric current vectors in the torsional Alfven wave of Figure 1.11 



3.26. Magnetic Reynolds numbers in a star 

When the dominant non-ideal effect is electrical resistance due to Coulomb interaction 
of the electrons with the ions and neutrals, an ionized plasma has a magnetic diffusivity 
of the order (the 'Spitzer value') 

with temperature T in Kelvin. A magnetic 'flux tube' of radius r = 10* cm and field 
strength 10^ G lies at the base of the solar convection zone, where the temperature is 
2 10^ K, and the density p = 0.2 g cm~'^. 

Calculate the magnetic Reynolds number for a convective flow of f ~ 100 m/s across 
the tube. Calculate the order of magnitude of the Hall drift velocity fn of the electrons 
in this tube (eq. |2.27 ), and from this a 'Hall- Reynolds number' Rh = v/vy{- Which 



direction does the Hall drift flow? 
Repeat the calculations for a magnetic filament of width 300 km and field strength 1000 
G in the penumbra of a sunspot, where the temperature is 5000 K, the density 10~^ 
g/cm~^ and velocities are of the order of 2 km/s. 
Conclusions? 



(Use Eq. 3.11 above only for estimates where a couple of orders of magnitude don't mat- 
ter. Quantitative values of the magnetic diffusivity require more detailed consideration 
of the physical conditions in the plasma.) 



3.27. Poynting flux in an Alfven wave 

Calculate the vertical component of the Poynting flux vector in the torsional wave of 
problem 3.24 What does the other component of S mean to you? 



3.28. Relativistic charge density 

Apply Lorentz transformations of the volumes occupied by the positive and negative 
charges moving in opposite directions due to a current density j. Show that in a frame of 
reference moving with uniform velocity v with respect to an electrically neutral plasma, a 



charge density appears that corresponds to (2.18). [If you want to make this easy, assume 



only the electrons are contributing to the current. To what extent is this justified?] 



3.29. Ambipolar drift 

A straight flux tube of circular cross section R is embedded in a field- free environment. 
Its field strength varies with distance r from its axis as B(r) = Bo[l — {r/R)"^]. Calculate 
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the current density in the tube. Show that the ambipolar drift velocity is irrotational. 
Which direction do the ions flow? 

In pressure equilibrium, the resulting pile-up leads to a pressure gradient between the 
charged component and the neutrals. In a stationary equilibrium, diffusion in this gra- 
dient just cancels the ambipolar drift. 
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4. Appendix 



4.1. Vector identities 

a-(b X c) = c-(a x b) = b-(c x a) (4.1) 

ax(b X c) = (a-c)b- (a-b)c (4.2) 

V x W = (4.3) 

V-(V X a) = (4.4) 

V X (V X a) = V(V • a) - V^a (4.5) 

V-(V'a) = • a + a- VV' (4.6) 

Vx(^a) = V'0xa + ■^V X a (4.7) 

V-(ax b) = b-(V X a) -a-(V X b) (4.8) 

Vx(a X b) = a(V • b) - b(V • a) + (b • V)a - (a • V)b (4.9) 

When memory fails these identities and others are rederived quickly using the properties 
of the permutation symbol or completely antisymmetric unit tensor: 

eijk, (4.10) 

where the indices each stand for one of the coordinates xi,X2,X3. Conventionally, in 
Cartesian coordinates, xi = x, X2 = y, x^ = z. It has the properties 

ei23 = 1, ejik = —eijk, ^jki = e^-fe, €ijk = if i = j. (4-11) 

It remains unchanged by a circular permutation of the indices. A product of two symbols 
with one common index has the property 

where Sij is the Kronecker delta, and the sum convention has been used, implying sum- 
mation over a repeated index: 

3 

aA = ^ aibi = a • b. (4.13) 

i=l 

The cross product and the curl operator can be written as 

{a X h)i = eijkttjbk, (V x b)j = eijfc — . (4.14) 
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4.2. Vector operators in cylindrical and spherical coordinates 

(See also appendix in Tritton 1992). Let a, b be vectors and / a scalar quantity. 
Cylindrical coordinates (r, cf), z) : 



I d , df^ 1 
V / = {r ) + + 

V dr dr dcjP' dz^ 


(4.15) 


V.a=--(ra,) + -^ + ^ 
T Br r B(h 3z 


(4.16) 


(V X a) - ^— - ^ 

^ r d(f> dz 


(4.17) 




(4.18) 


fV X a) ^ fra ) ^ ^'^^ 
^ r dr r dcf) 


(4.19) 


(V^a r = V^ar -^-^^ 
r^ r-^ 


(4.20) 


(V a)^ = V + ^2 


(4.21) 


(V^a)^ = V^a^ 


(4.22) 



where V^a^ denotes the result of the operator acting on regarded as a scalar. 

{ai-Vh)r = ai-Vhr-h^a^/r (4.23) 

(a-Vb)0 = a- V 60-6^0,. /r (4.24) 

(a-Vb)^ = a-V6;^ (4.25) 

where a-V 6^ denotes the result of the operator a-V acting on hi regarded as a scalar. 

Spherical coordinates {r,0,4>): 

V / = -,Q-^ir -) + ^^(-n^^) + (4.26) 

V-a = ~{r'ar) + ^^(a^sin^) + (4.27) 
r^ or r sm 6 ou r sm 6 o(p 

(V X a), = ^ (a^ sin 6)-^ (4.28) 
r sm t' at' ^ o4> 

N 1 5a,. 15,, 

V X a e = Trf - - TT ^"-^ ^"29 
r sm t) 0(p r or 

(Vxa), = l[|;(ra^)-|^] (4.30) 

(V^a), = V'ar - 4k + (a, sine) + (4.31) 

sm y oa sm Oip 

(V^a). = V^c.. + ^[^-^-^^l (4.32) 
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4.3. Useful numbers in cgs units 



Abbreviated here to 4 digits. (For recommended accurate values see 



|http: / /physics.nist.gov/cuu/Constants ) 



speed of light c = 2.998 10^° cm s'^ 

elementary charge e = 4.803 10^"'^'^ (erg cm)^/^ = g^/^ cm^/^ s~ 

gravitational constant G = 6.674 10^^ cm^ s~^ g"'^ 

proton mass rrip = 1.673 10^^'^ g 

electron mass rrie = 0.911 10^^^ g 

Boltzmann constant k = 1.381 10~^^ erg 

gas constant IZ (?» k/nip) = 8.314 10^ erg K^^g~^ 

Stefan-Boltzmann constant a = ac/A = 5.670 10~^ erg s^^ cm^^K"^ 

electron volt 1 eV= 1.602 lO"^*^ erg 

astronomical unit 1 AU = 1.496 10^^ cm 

solar mass Mq = 1.989 10^^ g 

solar radius Rq = 6.963 10^*^ cm 

1 year « tt • 10^ s 



4.4. MKSA and Gaussian units 

The choice of Gaussian units, where E and B have the same dimensions, is informed 
by the knowledge that they are components of a single quantity, the electromagnetic 
tensor. It is a natural choice in astrophysics, where quantities measured in earth-based 
engineering units are less likely to be needed. 

The MHD induction equation is the same in MKSA and in Gaussian quantities. In the 
equation of motion the only change is a different coefficient in the Lorentz force. In 
MKSA: 

Lorentz force: Fl = (V x B)xB//io 

Magnetic energy density: Cm = -B^/(2/io) 

The MKSA unit for the magnetic induction i?, the Tesla, is worth 10^ Gauss. The 
MKSA forms of E and j in MHD: 

Electric field strength: E = — (v x B) 
current density: j = (VxB)//io 

In Gaussian units the difference between magnetic field strength and induction, electric 
field and displacement, are dimensionless factors given by the relative permeability /ir 
and relative permittivity of the medium. In vacuum both are equal to unity, as 
assumed here. They have to be taken into account in some high-density environments, 
however. 
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